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Abstract. We construct explicitly (non-polynomial) eigenfunctions of the 
difference operators by Macdonald in case t = q k , k £ Z. This leads to a 
new, more elementary proof of several Macdonald conjectures, first proved by 
Cherednik. We also establish the algebraic integrability of Macdonald opera- 
tors at t = q k (k S Z), generalizing the result of Etingof and Styrkas. Our 
approach works uniformly for all root systems including BC n case and related 
Koornwinder polynomials. Moreover, we apply it for a certain deformation of 
A n root system where the previously known methods do not work. 



1. Introduction 

One of the goals of this paper is to present a new (essentially, "non-polynomial") 
approach to Macdonald polynomials. These polynomials were introduced in late 



80-s by I. G. Macdonald in [Ml] as, informally speaking, a "discrete spectrum" of cer- 
tain remarkable symmetric difference operators. Since these (Macdonald) operators 
are self-adjoint with respect to a specific scalar product, Macdonald's theory leads 
directly to the families of multivariable orthogonal polynomials. From that point 
of view, they generalize various classical orthogonal polynomials of one variable. In 
fact, there are as many families of Macdonald polynomials, as simple complex Lie 
algebras, or root systems. Each family depends, apart from a root system, on (at 
least) two parameters q, t and specializes to several remarkable families of symmet- 
ric functions. Among them are Schur functions and characters of corresponding 
Lie groups, Hall-Littlewood functions, Jack polynomials or, more generally, mul- 



tivariable Jacobi polynomials due to Heckman and Opdam [HO|. All this makes 
Macdonald polynomials very interesting from the point of view of the representa- 
tion theory, combinatorics, special function theory and mathematical physics. This 
also makes clear that Macdonald polynomials are highly non-trivial. Thus, it is not 
surprising that their various properties formulated by Macdonald as conjectures, 
remained unproven for quite a long time. A remarkable progress has been achieved 
by Cherednik, who proved the so-called norm and evaluation conjectures and the 



symmetry identity for all root systems [CI, C2[. Cherednik's approach is based on 



his theory of double affine Hecke algebras and it remains one of the major achieve- 
ments in this area. As an introduction into Cherednik's theory we recommend to 



the reader a nicely written survey by Kirillov [Kil . 

One of the results of the present paper is an independent proof of these three 
Macdonald's conjectures. Our approach uses some remarkable properties of the 
Macdonald operators in case t = q k with integer k. Ideologically, it goes back to the 
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paper by Veselov and the author [CV], where the quantum Calogero-Sutherland- 
Moser problem was considered for some special values of the coupling constant. 
Recall that the Calogero-Sutherland-Moser problem [Ca, Su] describes N particles 
x\, . . . , xn on the line whose pairwise interaction is given by the potential u(x) = 



m(m + 1) sin" 
operator in 1 



JV. 



In the quantum case its Hamiltonian is the following Schrodinger 



H = —A + m(m + 1) sin 2 (xi — xj) . 
i<j 

It is a celebrated example of a completely integrable TV-body problem, and there 
are quite a few exact results about it. In particular, it is completely integrable 
in Liouville sense, i.e. H is a member of a family of N commuting partial dif- 
ferential operators (quantum integrals) H± = H, H2, ■ • • , -ff/v- Moreover, as it was 
demonstrated in [|CV| , CSV |, for special values of the coupling constant m this 
problem becomes " much more" integrable. Namely, it turns out that the Calogero- 
Suthcrland-Moser problem for m G Z is algebraically integrable, i.e. its quantum 
integrals Hi, ... , Hn are a par t of some bigger commutative ring of partial 
differential operators (see [ CSV | for precise formulations and results). Moreover, 
Spec(9\) is an affine algebraic variety whose points parameterize Bloch eigenfunc- 
tions of H. This is a multidimensional analogue of a phenomenon, well known 
from the finite-gap theory in dimension one [DMN, Krl]. The main difference with 
the one-dimensional case is that the corresponding algebro-geometrical data are 
very "rigid" and really exceptional, which makes the existence problem for related 
multidimensional Baker- Akhiezer functions extremely difficult. This difficulty was 



overcome in [CSV] with the help of the theory of multivariable hypergeometric 
functions due to Heckman and Opdam [HO|. Recently, a direct independent proof 
was obtained in [Chl|. 

Continuing |Chl|, here we demonstrate that a similar phenomenon appears for 
H being replaced by any of the Macdonald difference operators, namely, they all 
are algebraically integrable for special integer values of the parameter(s). Note that 
in case of R being of A n type the Macdonald operators coincide (up to a certain 
gauge) with the trigonometric version of the elliptic Ruijsenaars operators from fRj, 
introduced as a generalization of the elliptic Caloger-Moser problem. We should 



mention at this point the paper by Etingof and Styrkas ES| where the algebraic 



integrability has been established for the Macdonald-Ruijsenaars operators with 
t = q k , k € Z. Their approach was based on earlier results by Etingof and Kirillov 



|EKI , EK2 ] who gave an interpretation of Macdonald polynomials for R — A n in 
terms of the representation theory for quant um groups. This delivers independent 
proofs of several results in case R — A n , see [EK3, EK4, Ki2 |. It is worth not icing 
that the symmetry identity in this case has been proved first by Koornwinder [ Kol| 
(see chapter VI o f Macdonald's book [M4Q. However, his method, as well as the 
methods of [EK2, ES|, does not extend to other root systems. 

One of the advantages of our approach is that it works equally well for all root 
systems and (we believe) is simpler comparing to [CI, C2[. The main object is what 
is natural to call a Baker- Akhiezer function ip for Macdonald operators. In case 
R = A n this coincides with the -^-function from JES[ . In that part which goes back 



to the papers [CV, CSV], our considerations have very much in common with [ES|. 
The main new ingredient is a fairly elementary construction of the i/>-function. Our 
main observation (Proposition 2.1 ) is that the Macdonald operators in case t = q k , 
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k € Z_ , act naturally in the coordinate ring of a certain very specific afHne algebraic 
variety. This implies the existence of ip which is our central result. From this we 
derive the duality, which reflects a certain symmetry between the two arguments x, z 
of ^-function. Basically, it means that, as a function of x, ip is an eigenfunction 
of the Macdonald operators related to the root system R , while, as a function 
of z, it is an eigenfunction of the Macdonald operators related to the dual root 
system. Thus, we observe on this level the bispectrality of Macdonald operators, if 
one uses the terminology of the fundamental work by Duistermaat and Griinbaum 
on bispectral problem [DG]. It is worth mentioning that in contrast with the one- 
dimensional case, bispectral problem in higher dimensions is much less investigated. 
For an interesting example related to Knizhnik-Zamolodchikov equation see recent 
paper |TVfl . 

Notice that our proof of the existence of ip is an effective one and gives a closed 
expression for it. The resulting formula generalizes the main result of ||Chl and it is 
a discrete version of one remarkable formula by Berest, who found in [Be] an elegant 
"universal" expression for the ^-function for the quantum (rational) Calogero- 
Moser problem. His derivation, however, was based on a crucial assumption that 
such a if) does exist. As we mentioned above, that type of existence problems 
is highly nontrivial in dimension > 1. Remarkably enough, our approach being 
inspired by the Berest's result, allows us to do these two things simultaneously: we 
prove the existence of i/j by a direct derivation of a discrete analogue of the Berest's 
formula. 

Everything extends without much difficulties to BC n case. A proper general- 
ization of Macdonald's theory for this case was suggested by Koornwinder [Kol|. 
The resulting orthogonal polynomials (Koornwinder polynomials) can be viewed 
as a multivariable analog of the celebrated Askey- Wilson polynomials [AW]. Van 
Diejen [vDl] showed that Koornwinder polynomials are joint eigenfunctions of n 
commuting difference operators for which he gave an explicit expression (initially, 
Koornwinder constructed one operator only). Further, in [vD2| the symmetry iden- 
tity (or self- duality) was established for a certain subfamily of Koornwinder poly- 
nomials. Then, finally, Sahi [^j proved duality in general case, using a proper 
generalization of Chercdnik's double affine Hecke algebra. Together with earlier 
van Diejen's results |vD2[ ] this implied also the evaluation identity and the norm 
identity, conjectured by Macdonald. Our approach leads to an independent proof 
of these results. As above, the key ingredient is the algebraic integrability of the 
Koornwinder operator for special integer values of the parameters. 

One of our primary motivations for this work was, in fact, our at tempt to find 
a difference version of the deformed Calogero-Moser problem from [CFV1]. It is 
related to what can be viewed as a one-parameter deformation of A n root system. 



In [ Ch2 , guided by duality, we were able to find a rational difference version of 
that quantum problem. Here we consider its trigonometric version, proving its 
(algebraic) integrability. Similar to the usual A n case, the constructed difference 
operator is self-dual: corresponding Baker- Akhiezer function ip(x, z) is invariant 
under permuting x, z. A natural elliptic version seems to be integrable, too. We 
would like to stress that while for the root systems the approach based on affine 
Hecke algebras seems to be the most adequate and powerful, in the deformed case 
none of the previously known methods can be applied (at least, straightforwardly). 
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Thus, it would be very interesting to find an algebraic structure which underlies 
that "deformed" root system. 

The paper is organized as follows. In Section 2 we recall the definitions of 
difference operators and polynomials due to Macdonald. Then we make our central 
observation about Macdonald operators in case t = q~ m (m £ Z+). The following 
two sections form the core of the paper. 

In Section 3, we define a Baker-Akhiezer function ip(x, z) associated to a datum 
which consists of a root system R and some additional integer parameter (s). %j) 
is determined by prescribing its analytic properties in the ^-variable, which is a 



proper modification of the approach from ]CV| E£]. We prove first its uniqueness 
up to a normalization and then the existence, obtaining as a by-product a discrete 
version of the Berest's formula. 

In Section 4, we explain how one should normalize tp to achieve a remarkable 
symmetry between x and z. This leads us directly to the duality theorem, which is 
the central result of this section. 

In Section 5 we derive various corollaries from achieved results. First, we obtain 
algebraic and complete integrability of the Macdonald operators and prove the 
existence of the so-called shift operators. Then we explain how our ip relates to 
Macdonald polynomials, this generalizes Weyl's character formula and is similar 
to a relation between symmetric and nonsymmetric Macdonald polynomials (see 
|M3| ) . As a corollary, we observe a nice "localization" property for Macdonald 
polynomials in case t — q~ m [m € We conclude the section by explaining how 
our results lead to a proof of the norm identity, evaluation formula and duality for 
Macdonald polynomials. 

In Section 6 we sketch how to extend our approach to BC n case. This is parallel 
to the previous sections, so we skip most of the proofs. The main difference comes 
from n = 1 case which is technically more difficult compared to Ai. 

Finally, in Section 7 we discuss a generalized Macdonald-Ruijsenaars model, 
related to the deformed A n system. 

Acknowledgments. I am grateful to Yu.Berest and A.P.Veselov for stimulating 
discussions. The work was supported by EPSRC (grant GR/M69548). 

2. Difference operators and polynomials by Macdonald 

2.1. Notations. Let V = K™ be a Euclidean space with the scalar product denoted 
as (it, v). Consider an arbitrary root system R € V which is, by definition, a finite 
set of vectors (roots) a € V with the following two properties: 
(1) Va S R the orthogonal reflection s a in V 



2- -a 

(a, a) 



leaves R invariant, s a (R) = R; 
(2)Va,(3eR 2^§eZ 



' (a, a) 

(see |B|] for the details). 

The second property implies that the so-called root lattice Q generated over Z 
by the roots a G R is invariant under all the reflections s a and, therefore, under 
the whole Weyl group W generated by s a , a G R. The vectors a v = 2a/ (a, a) 
form the dual root system i? v and we denote by Q v the lattice generated by all 

4 



a v E R v . Introduce also weight and coweight lattices P, P v as 

P = {tt e V\ (it, a v ) S Z Va v G R v } (weights) , 

P v = {tt £ V\ (tt, a)eZ Va € R} (coweights) . 

From the definitions one has inclusions Q C P, Q v C P v . Taking Z>o instead 
of Z in the last two definitions leads to dominant weights (coweights) P + and PY , 
respectively. 

Let us fix some basis of simple roots a\, . . . ,a n in R, this determines a decom- 
position of R and i? v into positive and negative parts: 

R = R+ U (-E+) , i? v = i?^ U (-RX) ■ 

The elements of the basis, dual to a( , . . . , o;„, (u>i,Oj) — Sij , are called the 
fundamental weights for the system i?. Similarly, one defines the fundamental 
coweights b^. (pi, aj) — 5^ . 

In these terms the root lattice Q is simply 



Q = Za 4 



Its positive part Q + is obtained by replacing Z by Z + . Similarly, the cone of the 
dominant weights P+ is 



P+=0Z + c 



»=l 

Below we will assume that i? = A„, . . . , Gi is reduced and irreducible. The case 
R = BC n is considered in section 6. 

2.2. Macdonald operators. Let R be a (reduced, irreducible) root system in 
Euclidean space V . Let us fix q 6 C x and a set fc of VK-invariant parameters 
fc Q = k a v , i.e. fc Q = k wa for any a £ i? and w £ W. Below we will sometimes use 
the related parameters t a = t a v := q ka , denoting by t = q k the whole set {i a }[] 
Throughout the paper we will suppose that q is not a root of unity. Something still 
can be done in case of roots of unity, though we will not touch these issues here 



(see (Ki2|, |Cg]). 

For v € V we denote by T v the operator acting on a function f(x) as a shift by 
v in x E V: 

(TV) (a;) = /(! + «). 

Later, we will deal with the functions depending on two variables x, z € V, in that 
case we will use subscripts to distinguish between shifts T£ and T", acting in x 
and z, respectively. To introduce Macdonald operators, we need the notion of a 
(quasi)minuscule coweight. 



Definition. (1) A coweight tt G P v is called minuscule if —1 < (tt, a) < 1 for all 
a e R. 

(2) A coweight tt S P v is called quasiminuscule if it belongs to i? v and — 1 < 
(it, a) < 1 for all a G i?.\7r v . 



lr To define g fc for arbitrary k we fix the value of logq, so g fc := e kl ° sq . 
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Using tables from Q one can check that all root systems except E$, F4, G2 have 
at least one nonzero minuscule coweight. Meanwhile, for any root system R the 
coroot n — 9 V where 9 is the maximal root for R, will be quasiminuscule (see Q| ) . 

Now let 7T e P v be a minuscule coweight for the system R. Then the cor- 
responding Macdonald operator is a difference operator in x € V defined as 



follows [Ml, M2 



+ n (a,x) _ ±—l„—(a,x) 

d *= e ^ tt < «^)= n g o( a,,)_ ;-(l„) ■ 



(o>,t)>0 



For a quasiminuscule coweight tt the formula is slightly more complicated: 

D"= «r(T T -l)+ J] ^ 2(P ' T) ' (2-2) 



r£Wn r£Wn 



where 



jwf */ TT '-a'f ''a *f /<j 0*1 

a (a,x) _ fl -(a,x) 11 n l+(a.x) _ -l-(a,x) ^•°>' 



(",t)>0 (a,r)=2 



and 



QG-R+ 



Remark. Note that in the formula (2.3) for a T the last product contains one factor 



only (with a — r v ). Written this way, it formally makes sense for a m inuscule 



coweight, too. Indeed, a formal substitution of a minuscule 7r into (2.2) leads to 
( pTl| ) because the constant term ^ q- 2 (p> T ) — ^ a T will be zero in that case. In the 
remainder of the paper, we will refer to the formula (|2.2|) since it covers both cases. 



Remark. The following function Afc(x) plays an important role in Macdonald's 
theory: 

00 1 _ n 2i+2(a.x) 

A = A k (x) = J] II ^T2M • (2-5) 

Using it, one can present the coefficients a T of the Macdonald operator D v as 
follows: 

a, = T W (A)/A, a^x) =a ff (^- 1 a;). (2.6) 



Example. In case i? = = {± (e* — ej) |i < j} C R™ with k a = k each funda- 

mental coweight 7r s = = e± + - ■ - + e s (s = 1, . . . , n) is minuscule and the corresponding 
operator D s — D^ 3 has the form 

D s = T TT" ._ g _ . V , (2.7) 

/c{i,...,n} ie/ y y 

|I|=S 3(1 

where T 1 stands for Y\ ie iT ei . These operators coincide (up to a certain gauge) 
with the quantum integrals of the trigonometric Ruijsenaars model introduced in 



2.3. Macdonald polynomials. The starting point for Macdonald's theory Ml 



M2 is that the operators , introduced above, preserve the space spanned by 



IF-invariant exponents, or orbitsums m x : 

m A (x) - ]T q 2 ^ , (2.8) 

rewx 

where A 6 P + is a dominant weight and the summation is taken over its VF-orbit. 
Moreover, the action of D" is lower-triangular: 

D*m x = cxxmx + ^2 c ^ m " ( 2 - 9 ) 

i/^A 

where the coefficients c Xv depend on q, t and v -< A means that A — t> belongs to 

Q+- 

To introduce the Macdonald polynomials, let us first agree about terminology. 
Throughout the paper by a polynomial in x we will always mean a function f(x) of 
the form 

From algebraic point of view this corresponds to considering a ring of Laurent 
polynomials in AQ = q^ uX \ where cj, are the fundamental weights. As well as a 
standard polynomial ring, it has a unique factorization property with q 2 ("' x *> being 
the only invertible elements. 

Now the Macdonald polynomials Pa = P\ (x; q,t) can be defined uniquely as 
polynomials of the form 

P\ = m A + a ^ m » > A e P + > ( 2 - 10 ) 

i^A 

which are eigenfunctions of D 7 *: 

D«P X = c xx P x . (2.11) 



The coefficients a Xv in ( 2.10 ) are rational in q, t and the polynomials Pa (A S P+) 
have a number of remarkable properties. In particular, they are orthogonal with 
respect to the following scalar product: 



(/, 5 ), = CT(/ 5 A fc A fc ) , (2.12) 



where CT means the constant term, A/~ is the function (2.5) and the bar acts on a 
function f(x) as f(x) = /(—as). This scalar product can be reinterpreted in terms 
of a certain integral, which makes perfect sense for noninteger k, too.m See |M1| 
for the details. 



Remark. It is not difficult to see that the eigenvalues caa for A e P + have the form: 

caa= £ q^+"\ (2.13) 



where p = pk is given by (2.4) 
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Remark. P\ is correctly defined if the diagonal terms in the action (2.9) are distinct: 



c\\ ^ c vv for all v<\. ( 2 -14) 
This is true for generic t a and in this case the polynomials ( 2. 1C ) are uniquely 



determined by (2.11 
For instance 



Their coefficients, however, have singularities at certain q, t. 
in case k a = —m a with m a G Z + which will be central in further 
considerations, some first P\(x;q,q~ m ) are not well-defined. Nevertheless, even in 
this case the condition (2.14) holds for sufficiently large A, i.e. if (A, a) > 2m a for 
all a > (at least, (2.14) will be true for a proper D w or their linear c ombi nation, 
sec | Ml , M2 1 ) . This means that for such A the w hole set of equations ( [2.11 ) for all 
(quasi)minuscule coweights together with ( 2.10 ) determines Pa correctly. 



2.4. Macdonald operators in case k a G Z_. Let us concentrate now on the case 
of integer multiplicities k a . It is known that some results in Macdonald theory are 
easier to achieve for integer k a , extending then to non-integer values by a proper 
"analytic continuation" in k a . However, instead of a common assumption k a G Z+, 
we will rather consider the case k a — — m a where m a G Z + . As we mentioned 
above, in this situation the Macdonald's theory is not complete. Nevertheless, 
the structure of the eigenfunctions of the operators can be described quite 
effectively Next proposition will be crucial for us. 

For given root system R and integer multiplicities m = {m a } introduce a ring 
91 which consists of all polynomials f(x) with the following properties: for each 
a G i?+ and j = 1, . . . , m a 



I ( * + \^ j - / 



1 • x 

x la 

2 J 



for q 2 ^ = 1 



(2.15) 



Proposition 2.1. Let D 7 * be a Macdonald operator , defined in accordance with the 
formulas (2.1)— (2.3). Suppose that all t a have the form t a = q~ m " with m a G Z+. 



Then the operator D n preserves the ring (2.15): D^^PC) C 91 



To prove this, we will look first at the rank-one case, R C V — K 1 . In this case 
we will denote by T the unit shift: (Tf)(x) = f(x + 1). Similarly, T s will stand for 
the shift by a scalar s. Let us consider a difference operator L of the form 

L = a{x)T + b{x)T- 1 . (2.16) 

Suppose that its coefficients a, b are meromorphic with simple poles at x = and 
with no other poles at x G Z. Further, let us fix an integer m G Z + and impose the 
following conditions on a, b: 

res x=0 {a + b) = 0, (2.17) 

a(j)=b(-j) fori =±l,...,±m, (2.18) 

o(m) = 0. (2.19) 

Introduce a ring 5Ho which consists of all meromorphic functions f(x) with no poles 
at x G Z and such that 

f(j) = f(-j) forallj = l,...,m. (2.20) 
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Lemma 2.2. Under conditions (2.17)— (2.1£ ), the operator L preserves the ring 
( pC| ): 

L{m a )cm a . 



Proof. First of all, for any / € £Ho its image Lf will be nonsingular at i e Z. 
Indeed, the only apparent pole is x — 0. However, it disappears since the residues 
of a and b are opposite and (T — T~ 1 )f\ x= o — due to ( pc|) at j = 1. 



Now let us prove that L/ still satisfies the conditions ( [2.20 ), i.e. that (T^ — 
T~i)Lf is zero at x — 0. A simple calculation gives us that 

( T i-T-j)Lf\ x = = a(j)f(j + l)-b(-j)f(-j-l) + b(j)f(j-l)-a(-j)f(-j+l) 
= a(j) (f(j + 1) - f(-j - 1)) + b(j) (f(j - 1) - f(—j + 1)) 



(here we used the conditions (2. IS)). The last expression must be zero for all 
j = 1, ... ,m due to @.18| ) and fl2.19| ). □ 



Our next lemma is a modification of the previous one for a three-term difference 
operator 

L = a{x) (T 2 - 1) + b(x) (T- 2 - 1) . (2.21) 

Now a and b may have simple poles at x = 0, —1 and at a; = 0, 1, respectively, 
with no other poles allowed at x G Z. Further, we fix an integer m as before and 
impose the following conditions on a, b (in case m > 1): 

resz^a + b) = 0, res x =-i(a) + res x= i(fe) = , 
a(j) = b(-j) for j = l,±2,...,±m, 
a(m) = a(m — 1) = . 

In case m = 1 these conditions must be replaced by the following: 

res x=0 (a) = res x=0 {b) = , res x= _i(a) + res a!= i(i») = , 
o(l) = 6(-l) = 0. 



Lem ma 2.3. Under the assumptions above, the operator ( 2.21 ) preserves the ring 
( PC| ): £(Wo) C • 



Proof. Take any / £ 5R , then the only possible poles of I// are x = ±1 and x = 
(if m > 1). We have : 

re S:c= o(£/) = re S:c= o(a) [/(2) - /(0)] + re S;c= o(&) [/(-2) - /(0)] . 

This is zero for m > 1 because /(2) = /(— 2) and res x= o(a + 6) = 0. Further, 

res :c= _ 1 (L/) - res^.^a) [/(l) - /(-l)] = , 
res x= i(L/) = re S;E=1 (6) [/(-l) - /(l)] = . 

So, L/ has no singularities at x £ Z. 



Now let us check that Lf still belongs to <K , i-e. (T j - T~*) Lf\ x=Q = for 
j = 1, . . . , m. A straightforward calculation gives us the following: 

(T j - T~ j ) Lf 

= a(x + j) [f(x + j + 2) — f{x + j)] + b(x + j) [f(x + j - 2) - f(x + j)] 
- a(x - j) [f{x -j + 2)- fix - j)} - bix - j) [fix - j - 2) - fix - j)} . 
For j > 1 each term in this expression can be evaluated at x — 0: 
( T i _ r -i) Lf \ x=o = a{j) [f{j +2 ) - /(j)] + [/(j - 2) - /(j)] 
- o(-i) [/(— j + 2) - fi-j)] b(-j) [fi-j -2)- fi-j)] 
= aij) [fij + 2)- fi-j -2)- fij) + fi-j)} 

+ bif) [fij -2)- fi-j + 2)- fij) + fi-j)} . 

Here we used the conditions a(j) = b(—j) and a{—j) — b(J). The resulting ex- 
pression for all j = 1, . . . , m will be zero due to the properties ( 2.20 ) of / and the 
condition a(m) = a(m — 1) = 0. 
Finally, for j = 1 we have: 

(T - T- 1 ) Lf = aix + 1) [fix + 3) - /(z + 1)] 

- bix - 1) [fix - 3) - /(z - 1)] 

+ iaix - 1) + 6(1 + 1)) [/(z - 1) - fix + 1)] . 

The last term is zero at x = because res^o {a{x — 1) + 6(z + 1)) is zero and 
/(l) = /( — 1). The first two terms after evaluating at x — give 

a(l)[/(3)-/(l)]-6(-l)[/(-3)-/(-l)] . 

This is zero for m = 1, 2 since a(l) = 6(— 1) = in this case, and for m > 2 since 
a(l) = 6(-l) and /(3) - /(-3) = /(l) - /(-l) = 0. □ 

The following two lemmas are a direct corollary of the previous two. 

Lemma 2.4. Suppose that a difference operator L of the form ( |2.16| ) is invariant 
under the change of variable x — > —x, i.e. 6(x) = a(— x). Further, let us suppose 
that 

(1) a, b have simple poles at x = and no other poles at x G Z , 

(2) a(m) =0. 
Then L(SHo) C <K - 

Lemma 2.5. Suppose that a difference operator L of the form ( |2.2l| ) is invariant 
under the change of variable x — > — x, i.e. 6(x) = a(— x). Further, let us suppose 
that 

(1) a, 6 have simple poles at x = 0,-1 and x = 0, 1, respectively, and no other 
poles at x e 7L , 

(2) a(m) = a(m - 1) = 0. 

(In case m — 1 we replace it by the condition a(l) = but now allow poles at 
x = ±1 only.) 
Then L(9to) C <K - 



One can formulate the following inversions of lemmas 2.2 and 
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Lemma 2.6. Suppose we are in a situation described in lemma 2.2. Moreover, let 
us impose extra conditions on the coefficients a, b of the operator ( [2.16| ) as follows: 

res x=0 (a) 7^ 0, 

a(j) for j = 1, . . . , m — 1 (only in case m > 1) . 
Let / be an analytic eigenfunction for L: Lf = Xf. Then / belongs to the ring D\ . 



Lemma 2.7. Suppose we are in a situation described in lemma 2.3. Let us impose 



extra conditions on the coefficients a, b of the operator ( 2.21 ) as follows: 
res x= o(a) 7^ (only in case m > 1) , 
resx=-i(a) ^ 0, 

a (j) 7^ f° r J — 1, . . . , — 2 (only in case m > 2) . 
Then each analytic eigenfunction / of the operator L must belong to the ring £Ho- 

Proof. Both lemmas can be proven by reversing the arguments used to prove lem- 



mas 2.2, 2.3. Indeed, since / is non-singular, Lf = Xf must be non-singular, 
too. Looking at the residues of Lf, we obtain the condition /(l) = /(— 1). Thus, 
Lf = Xf must also satisfy this condition, which gives more conditions on /, and so 
on. 

□ 



Let us apply all this to the difference operators 



D 1 = q - Q - T + ^ Q - T- 1 (2.22) 

„.T ~ — X ~X ~ — X v ' 



and 



2 (q x - q- x ){q x + x - q- x ~ x ) 1 ' 

i — x— -m j+mU„- x— m+1 rl x+m—l\ 

+ - rr-^ ^ >-(T- 2 -l). (2.23) 

which are the Macdonald operators in cases R = A\ = {±1} C K 1 and n = 1 
and 2, respectively. It is obvious that D\ and Z?2 satisfy the conditions of lemma 
2.4 and ^5 (provided that q is not a root of unity). Hence, they preserve the 



properties (2.20). Moreover, instead of x — we may consider any point x = xq 
with q x ° — q~ x ° = 0. Indeed, for such xo the symmetry x — > 2xo — x still does not 
change the operators D\ and L>2 (this reflects their invariance with respect to the 
corresponding affine Weyl group). So, after shifting the origin to xq one gets the 



operators with the same properties as in lemma 2.4 and ^5|. Thus, the operators 
Di, £>2 preserve, in fact, a bigger ring, which is an affine version of the ring (2.20). 

Namely, let us consider the ring of all analytic functions f(x) with the following 
properties: 

fi x + j) = f( x - j) for each j = l,...,m and q 2x — 1 . (2.24) 
The arguments above prove the following result. 



li 



Proposition 2.8. The operators (2.22), (2.23) with m S Z + preserve the ring 91 
of analytic functions with properties (2.24). 



This is, essentially, the rank-one case of Proposition 2.1. Now we are ready to prove 
it in full generality. 



Proof of Proposition \2.\ We will only consider the operator D n given by the for- 
mulas (|2.2j)-(^3[), since this covers the case (2.1), see Remark]!]. Choose any a € R, 
then we should prove that D n preserves the ring £H Q of functions with the following 
properties: 

f(x+\j<A =f(x-\jaA for j = l,...,m a and q 2 ^ = 1 . (2.25) 



Let s = s a € W be the reflection with respect to a. Since the orbit Wtt of the 
coweight 7r is s-invariant, it splits into pairs r, t' with r' = s(r) plus a number of 
s-invariant r's. This defines a splitting of D n into a sum of difference operators of 
one of the following three types: 

D = a r (x)(T r -l), s(t)=t, (2.26) 
Di = a T (x){T T - 1) +a T ,(x)(T T ' - 1) , t' = t - a v , (2.27) 
D 2 = a a v (T qV - 1) + a_ Q v (T~ qV - 1) ( only for quasiminuscule tt ) . (2.28) 

This follows directly from the fact that ir is (quasi)minuscule. Moreover, since D" 
was obviously VK-invariant, each of Di will be invariant under the reflection s. In 
particular, a T (x) = a T (sx) in Do, a T i(x) = a T (sx) in D±, and a_ Q v (x) = a a v[sx) 
in Z? 2 . 

We claim that each of Di preserves the ring ( 2.25 ). First, note that in case of 
Dq a T has no pole at q 2 ( a > x *) = 1 since it is s-invariant. So, a T itself belongs to the 
ring D\ a . Also it is clear that operator T T — 1 preserves this ring (the shift acts in 
direction, orthogonal to a). Hence, Do(D\ a ) C £H Q . 

Now let us consider D\, it has the form 

a T (x)T T + a T '{x)T T — a T (x) — a T >{x) . 

The sum a T + a T ' is s-invariant, hence, it is nonsingular at q 2 ^^ = 1 and belongs 
to y\ a . So, we may ignore it and consider the first two terms only. Further, note 
that we can present r, r' as 

1 V / V 

r = —a + t> , r = cr + v 

2 ' 2 

for a certain u such that (a,u) = 0. Hence, T 11 will preserve the ring 9\ a and we 
can also ignore it, reducing D\ to 

a r T3 QV +a T /T"2 aV . (2.29) 

This operator is s till s-invariant. Moreover, since (a, t) = 1 in this case, we see 
from the formula (2.3) for a T that 

a T (x + ^m a a v ) = for q 2{a ^ = 1 . 

Now in absolutely the same way as it was for the operator ( 2.22p , we conclude that 
( 2.2SD preserves the ring 

In the same manner the case of D 2 reduces to (2.23). □ 
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Remark. One can show that the ring is finitely generated, therefore it can be 
viewed as the coordinate ring of a certain quite specific affine algebraic variety. For 
instance, for A± case it is a rational curve with m double points. 



Remark. In our proof of Proposition |2 . l| essential ingredients were the VF-invariance 
of the operator D" and specific location of poles and zeros of its coefficient s. Thi s 
has ce rtain parallels with the residue co nstruc tion of Hecke algebras from [ GKV |, 
BEG]. Moreover, using the results from | BEG |, one can prove that for k a = —m a g 
Z_ all Macdonald- Cherednik operators Z) 1; . . . , D n coming from VF-invariant part 
of the double affine Hecke algebra |C1| will preserve the ring 9L All these operators 
will commute with D n . However, as we will see later, the centralizer of D n in case 
k a G Z_ is much bigger, and it contains many non-symmetric difference operators, 
preserving £H. 



One can formulate an analogue of Proposition 2.1 for the case of positive k a G 
Z + , too. This is because these two cases are related through a simple gauge trans- 
formation. Let 7r be a (quasi)minuscule coweight and D m denote the corresponding 
Macdonald operator (2.3) with t = q~ m . For m a G Z+ introduce a function 8 m (x) 
as follows: 



s m {x)= Yl Yl [(«.*) 



(2.30) 



aS-R+ j = — m a 



where [a] denotes [a] = q a — q 
checked by a direct calculation. 



The following fact is well-known and can be 



Lemma 2.9. For m a G Z + we have the following relation between the Macdonald 
operators D m and D m i with m! = — 1 — m: 



oD. n 



Corollary 2.10. LetSK = SH m be the ring (2.15). Then the Macdonald operator D, 
given by ^^-(^) with t = q m+1 , preserves the ^-module U = S^IR: D(U) C 
U. 



3. Baker-Akhiezer functions for Macdonald operators 

We keep mostly the notations of the previous section. So, we consider an arbi- 
trary (reduced, irreducible) root system R and fix a VF-invariant set m of multi- 
plicities m a G Z + . Ou r pu rpose is to construct eigenfunctions of the corresponding 
Ma cdonald operators (2^) with t a = g~ m ° (so k a — —m a in notations of section 
2.2) . Keeping this in mind, we introduce p = p{m) instead of p = pk = —p m from 



(2A) 



= p m = r ^ m a t 



ueR+ 

We will also use its counterpart for the dual root system i? v 



(3-1) 



(3.2) 



( Warning: p v ^ 2p/{p, p) !) 
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In this section we often will deal with functions of two variables x, z € V. We 
will keep calling a sum f(x) = X^ep q 2 ^' x ^ a polynomial in x. However, switching 
to z, we will also switch from the root system R to its dual i? v . For instance, 
by a polynomial in z we will mean a sum f(z) = 5Z„ g pv q 2 ^ v ' z h Such (perhaps 
confusing) terminology is caused by our implicit identification of the vector space 
V and its dual V*. To distinguish between these spaces, it would be natural to 
assume that R C V* and i? v C V. In this case x and z would lie in V and V* , 
respectively. However, we prefer not to do this, keeping things simple. 

We will apply the term quasipolynomial (in x or in z) to a function of the form 
q 2 ( x ' z ' f, where / is polynomial in x or in z, respectively. 

For a polynomial f{x) = X^eP o>vq 2 ^"' x ' by its support supp(/) we will always 
mean the convex hull of all points v with a u ^= 0. In the same way we define 
the support of g(z) — X^ppv a v q 2 ^' z \ For a quasipolynomial in x of the form 
<j> = q 2 ( x ' z > f(x) by its support supp(</>) we will simply mean the support of f(x). 
In case when <fr is quasipolinomial in x and z at the same time, it usually will be 
clear which support we are considering (either in x or in z). 

3.1. Baker Akhiezer function: definition and uniqueness. Let i/j(x,z) be a 

function of two variables i,z6 V of the form 

^ = q 2{x ' z) £ ^q 2 ^ , (3.3) 

where the coefficients i p u = ip v (x) depend on x, (x,z) is the scalar product in V 
and the summation in (3.3) is taken over all coweights v G P v from the following 
polytope TV: 

N = {v = p 1 - \ ^<l a < m a } . (3.4) 

Using our conventions about terminology, these conditions on a function ip can be 
rephrased as follows: "0 is quasipolynomial in z with supp(^) C M . 

Suppose that ip satisfies also the following conditions: for each a £ R and 
s = 1, . ..,m a 



ij) (x, z + 2 SC ^j = ^ (^ x > z ~ 2 Sa ) = ^ 



(3.5) 



Definition. A function ip(x,z) with the properties (3.3)-(p.q) is called a Baker 



Akhiezer (BA) function associated to the data {i?, m}. 

Our terminology is motivated by the fact that in case R = A\ such a ip is a 



Krichever's Baker-Akhiezer function [Krl, Kr2] associated to a specific singular 



rational curve. In contrast with the one-dimensional case, in higher dimension the 
main problem is to prove the existence of such a function. We do this in the next 
subsection. Let us presume now that such a ip does exist. 



Proposition 3.1. Properties (3.3)- (3.5) determine ip uniquely up to a factor 
depending on x. 

Proof is based on the following two lemmas. 
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Lemma 3.2. Let a quasipolynomial in z tp(x, z) = q 2 ( x ' z ) J2^eP v Vv<? 2( ' ly ' z ' ) satisfy 
the conditions fl3.5| ) . Then for each a v G R v and for any v E P v the set of integers 
j such that Vv+ja v 7^ either is empty or contains at least two integers j\, j% with 



Lemma 3.3. Let l\,...,l r be a set of non-parallel segments in affine Euclidean 
space V and SI c V be a convex domain in V. Suppose that for each li and for any 
line /, which intersects il and is parallel to k, the intersection SI n I has the length 
greater or equal than Then SI can be presented as S7 = Cl'#Af for some convex 
domain SI', where AT — lijfchjf 1 ■ ■ ■ jpl T - Here # denotes the Minkowski addition in 
V. 



We recall that the Minkowski sum of two subsets A, B of an affine space V is 
formed by all the points a + 6, where a and b run over A and respectively. The 
addition of points, of course, depends on the choice of origin, but the resulting set 
will be the same up to a shift. This operation is relevant to the multiplication 
of polynomials: if /i,/2 are two polynomials in x, and Ni — supp(/i) then N = 
supp(/i/2) is the Minkowski sum of N\ and N2. 



Proof of Lemma \3.q . For a given a € i?+ , substitution of the ip into 
following set of relations: 

J2 ^Q 2M (q {sa ' x+v) - q -( sa > x +^ = for q 2{a ^ z) = 1 
(s = 1, . . . ,m a ). 

These relations split up into separate linear equations for each "a-string" Vj — 

E^'(fe) S -(%r s ) =° (s=l...m a ), (3.6) 

3 

where ijjj := ip Uo +j a v and qj :— q 2 - 7 q( a > x +"o) . 

Suppose now that among the coefficients ipj only ipi,... ,ipm a do not vanish. 
In this situation we would have a homogeneous linear system of m a equations for 
m a unknowns ipj. Thus, it would be sufficient to show that this system is non- 
degenerate for generic x. To see this, we can look at the asymptotic behaviour of 
its determinant at large x (cf. |ES| ). More precisely, we consider the corresponding 

matrix A = (0^)1^=1 ma with aij — (qj) 1 — (qj)~ l where qj :— g 2j q( a > x + u o) _ Then 

for large x such that qi a > x + v o) ^> 1 the determinant det A asymptotically equals the 
Vandermonde determinant det(gj) which is nonzero since q is not a root of unity. 

Thus, the system is non-degenerate for generic x and all ipj must vanish, which 
proves the lemma. □ 



3.5) gives the 



Proof of Lemma \3.j[ . We will prove the lemma by induction in the number s of the 
segments. For an easier reference to the assumptions from the lemma, let us say in 
such a situation that J! dominates over the segments l\, . . . , l s . 

First, suppose we have just one segment l\ and a convex domain SI which domi- 
nates over l\. Then SI = Zi#Sl', where ST is the intersection SI n T(fl) of fl and its 
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image under the shift T for a vector l\ , associated with the segment l\ (in cither 
of two possible directions). This proves the lemma in case s = 1. 

Now suppose that £1 dominates over l\, . . . ,l s . Take the first segment l\ and 
consider the convex domain il' constructed above, so we have f2 = h#Q,'. We 
claim that f2' still dominates over h, . ■ ■ , l s . 

To prove this, let us take, for instance, I2 and choose any line parallel to it. 
Now we may consider a two-dimensional section of the f2 passing through / and 
parallel to 1%. The resulting two-dimensional domain will, obviously, dominate over 
li and Z 2 . Thus, essentially, we need to check the lemma in dimension two, for 
s = 2. This is very simple, and we leave it to the reader. 

So, we have proved that f2' dominates over I2, ■ ■ ■ ,l s - Now the statement of the 
lemma follows by an obvious induction. □ 



Proof of Proposition \3. 4 First, notice that the polytope JV in (3A) is exactly the 
Minkowski sum of the segments m a a v with a v £ i?Y (abusing notation, we denote 
a vector and associated segment by the same symbol). Now let ip be any BA 
function. By definition, we have an inclusion 

supp(V') C TV . 



On the other hand, Lemma 3.2 implies that the polytope supp("0) must dominate 
over each of the segments m a a y , a v G i?^. Hence, by Lemma |3.3| it must contain 
(a copy of) the polytope ( |3.4[ ) which is their Minkowski sum. Altogether this proves 
that for each (nonzero) BA function ijj one has the equality 

supp (■(/>) = N . 

Let now t/>', ip" be two Baker- Akhiezer functions. Consid er th ei r li near com- 
bination ip = ip' — c(x)ip" , which still satisfies the conditions (3.3 )-(3.5). We can 



choose c(x) in such way that the resulting function ip will have zero coefficient at 
one of the vertices of the supp(-0') = supp(-0") = N ' ■ So, we will have a strict 
inclusion supp^) C M. Thus, the only possibility is that such ip is zero, hence, 
ip' = c(x)tp" . □ 



Corollary 3.4. For a BA function ip, the nonzero coefficients ip u in (3.3) can ap- 
pear only for v = p v — X^agfl, laCt v with integer l a . In other words, the summation 
in (|3.3|) is taken effectively only over the set p v + Q v C P v . 



Proof. Suppose we have other terms, then let us remove them from the sum (3.3). 
This would not affect the conditions ( |3.5| ). Indeed, in the process of proving lemma 
3.2 we saw that these conditions split into separate linear equations involving j/'s 
from the same coset in P v /Q v . But from the uniqueness of ip it follows that the 
resulting function must remain the same. Hence, there were no other terms at 
all. □ 

16 



3.2. Existence of BA function. Comparison of Proposition |2.1| and conditions 



(3.5) suggests the idea to use a Macdonald operator acting in the z-variable in order 
to construct a Baker-Akhiezer function ip. Let lj G P be a (quasi) minuscule weight 
for the system R and be the Macdonald operator (2.2) corresponding to the 
dual system R v and acting in the z-variable: 

2(pV) 



E a T (Tl-l)+ J2 «" 



where p v is given by (|3.2[) and 



a T (z) 



n 

oE»: 
(a,T)>Q 



[(q v , z) - m a ] 



n 



a£»: 
(aV, T )=a 



[(q v ,z)-m a + l] 
[(a\z) + l] 



(3.7) 



(3.8) 



where [a] = q a — q~ a . 

Introduce the ring which is a counterpart of the ring (2.15) and consists of all 
polynomials f{z) = X^epv fyq 2tyV ' z ^ with the following properties: for each a G R+ 
and j = 1, . . . , m a 



for g 2 (° v .*) = l. 



(3.9) 



Then, according to Proposition 2.1, the operator will preserve the ring 91 v : 
D^(5K V ) C $H V . 

Now we need one technical lemma which shows that the action of D% on 5K V is 
" lower-triangular" . 



Lemma 3.5. Let D = be the Macdonald operator (3.7). Suppose that both / 
and / = Df are polynomials in z: 

m = E ^ 2{v,z) > fo) = E f^ 2(u ' z) ■ ( 3 - 10 ) 

v€P v i/€P v 

Then supp(/) C supp(/). Further, let A be a vertex of the polytope N = supp(/), 



then the ratio c\ = fx/ fx of the corresponding coefficients in ( 3.1C ) can be calcu- 
lated as follows. First, choose generic v S V such that (v, A) > (v, v) for all v G N 
and put 

1 



pi 



2 E 

agfl: (a,v)>0 



Then one has c A = £ TeH ^ 9 



2(r,A-pV) 



This can be proven similar to the proof of ( |2.9| ) and ( |2.13 ) in [Ml, M2|. A key 
point is that inclusion A#C C £?#C for convex domains A, i?, C implies A C B. 

A similar result is true for quasipolynomials. Recall that for a quasipolynomial 
4> = q 2 ( x,z 'f(z) its support, by our conventions, coincides with the support of /. 
Thus, we have the following analog of the lemma above. 



Lemma 3.6. Suppose that both <fi and 
q 2 



(z) = q 2 ^ Y <f> u q 2 ^ 



Dtp are quasipolynomials in z: 

q 2[x ' z) y M 2(l/ ' Z) ■ (3.11) 



v£P v 



17 



Then supp(0) C supp(</>). Further, let A be a vertex of the polytope N 
then the ratio c\ 
lated as 

2{r,x+\-pl) 



supp(^), 



4>x/^>\ of the corresponding coefficients in ( 3.11 ) can be calcu- 

(3-12) 



E 

t£Woj 



where p\ is defined in the lemma above. 



Now everything is ready to construct a BA function. The idea is very simple. 



We start from the quasipolynomial 4> = q 2 ^- x '^Q(z), where 



Q(z[ 



(3.13) 



where [a], as usual, denotes q a — q~ a . This polynomial is especially chosen to guar- 
antee that 4> satisfies the conditions (|j|) in z. Thus, applying D = successively 
to 4>, we will always get a quasipolynomial in z which will still satisfy the conditions 
( pTgj ) in z. Let us apply at each step an operator D — Ci, so 4>i+i = {D — Ci)4>i,4>o = 4>. 
The coefficients Cj will be adjusted to reduce supp(0i), see below. Finally, tp will be 
obtained after repeating this sufficiently many times. Before pro ceedi ng with more 
details, let us make one more remark. Notice that the formula ( 3.13 ) implies that 
the only nonzero terms in the initial quasipolynomial = q 2 ^ x ^ X^gpv 4> v q 2 ^ v ' z ^ 
are those with v e p v + Q v C P v . To see this, one should rewrite each factor 
[(a v ,z) + j] in fl3.f3| ) as q^ ' z \q j - q-^ 2 ^ »). The same will be true for all 
successive functions ipi. One can see this directly after rewriting in a similar way 
the coefficients (3.8) of the difference operator D. 



Now let us make everything more concrete. At the beginning we have 



A^o = supp(0 o ) = {v = p" 



laC 



m Q <l a < m a } . 



Further, as we mentioned already, each of 4>i will satisfy the conditions ( |3 . 5[ ) . Hence, 
due to lemmas 3.2 and 3.3, its support Nj C No must be a union of several copies 
of the polytope Af given by (3^). Let us fix generic v lying inside the positive 
Weyl chamber C. The linear functional (v, •) determines the height function on V. 
We will call a vertex of a polytope the highest (respectively, the lowest) vertex if 
it has maximal (respectively, minimal) height among all vertices. At each step we 
will choose the highest vertex A of the polytope Ni = supp(0;). For brevity, let 
us call the corresponding coefficient <p\ the highest coefficient of <p. Note that the 
initial polytope N Q (as well as all successive Ni) will be composed of the images of 
A^ under some of the shifts by v — X^agi?. '» qV with l a = 0, . . . , m a . If we look at 
the smaller polytope A/", then its highest (resp. lowest) vertex will be A = p v (resp. 
A = — p v ). Hence, the highest vertex of the polytope Ni must be of the form 



A = P V 



aeR+ 



la = 0, 



(3.14) 



Now we can kill the highest coefficient <p\ by applying D — ct where Ci = c\ is given 
by (3.12). Note that the vector p\ in the formula (3.12) will be simply p v (because 
A is the highest vertex). Thus, = c\ will be Cj = ^2 TeWuJ q 2 ( T ' x+l, \ or simply 

Ci = m w (x + v) 
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in terms of the orbitsum m w , 

On the other hand, let us look now what is happening with the lowest coefficient 
0_ p v . According to Lemma |3.6| , after application of D — a it gets the factor 
c_ p v — Cj. Note that for A = — p v the vector p\ in the formula (3.12) will be simply 
-p v . So, the formula ( |3.12[ ) gives us: 

c-pv = ^ g a (^)=tn & ,(a:). 

Obviously, c_ p v — Cj = m w (x) — m w (a; + v) is nonzero as soon as v ^ (recall that 
we assume that q is not a root of unity) . 

These considerations imply the existence of BA functions for all root system, 
which is one of our main results. 



Theorem 3.7. Let D — denote the Macdonald operator (3/7). Define ip(x,z) 
as follows: 



i> = Yl{D - nUx + v)) 



Q{z) 



(3.16) 



in accordance with the formulas ( p. 13] ), (3.15), where the product is taken over all 
v 7^ having the form v — ^2 aeR+ l a ct v with l a = 0, . . . , m a . Then 

(i) ip has the form (3.3)~(3.4); 

(ii) the coefficient ip- p v in its expansion ( |3.3| ) equals 



1>- f 



Y[(xa u (x) - m^ix + v)) ^ 0; 



(3.17) 



(Hi) ip is a Baker-Akhiezer function for the system R with multiplicities m — {m a }; 
(iv) as a function of z, ip is an eigenfunction of the Macdonald operator D: 

Dip = muj(x)ip . 



Proof. As we explained above, the constructed function ( [3.16] ) will be a quasipoly- 
nomial in z satisfying the conditions (3.5). The part (ii) follows immediately from 
the construction of ip. It implies that the polytope supp(V>) contains the polytope 
(3.4). On the other hand, the arguments above show that the highest vertex of 
suppV> must be A = p v , because every "higher" term ip p v +v has been killed after 
applying D — m iAj {x + v). Altogether, this gives us that supp(V') coincides with the 
polytope M in (3.4). Thus, the part (i) is also proven. 

Part (iii) of the theorem follows from the previous two and the remark above that 
ij) satisfies conditions ( |3.5| ). Finally, tp = (D — m UJ (x))^ must be quasipolynomial 
in z with all the properties (|3.3])-(3.5). Hence, it must be proportional to ip due 
to Proposition 3.1. However, application of D — m^^x) kills the highest coefficient 
ip p v, so ip must be zero, which proves the last part. □ 



Remark. Formula ( 3.16| ) is a trigonometric version of a related formula from | Chl |, 
which, in its turn, is a discrete version of the Berest's formula [Be]. 
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4. BlSPECTRAL DUALITY 



In this section we will explain how one should normalize BA function to achieve 
a certain symmetry between x and z variables. We start by looking closely at the 
rank one case. 

4.1. Ai case. For the rank-one case the existence of a BA function is a very simple 
fact, since in this case the number of "free" parameters equals the number of con- 



ditions (3.5). Let us consider the root system R = {±2} C K with m a = to S Z + . 
It is convenient to fix a scalar product on V = R as (u, v) = In this case we 

will have R = R v , Q = Q v = 2Z and P = P v = Z. We will fix R + = {2}, so 
p = p v = m. In accordance with (|3.3|), the Baker- Akhiezer function ip depends on 



two scalar variables x, z and has the following form: 



i> = <r Y\ (4.1) 



v— — m 



By definition, it must satisfy the following conditions: 

(T'i ~ T~ j ) i(j = for each j = 1, . . . , to and q 2z = 1 . (4.2) 

Similarly to (|3.6|), these conditions lead to the following linear system for the coef- 
ficients %b v : 

m 

J~] a,ijip- m+2 j = , i = l...m, (4.3) 

3=0 

where ay = (qj) 1 — {qj)~ % with qj := q- m + 2 i+ x . Introduce to x (to + 1) matrix 
A = (ciij), then the linear system above takes the form Av = 0, where v is the 
column v = (i/)- m , ip-m+2, ■ ■ ■ > V'm)- We know already that for generic x the matrix 
A has the maximal rank (equal to to), hence, its kernel is onc-dimcnsional. Using 
Cramer's rule, we find the values of ip v (up to a common factor): 

i>- m+ 2s= (-l) J 'detA«, 

where A^ s ' is obtained from A by deleting its s-th column. This gives us the ratio 
ib- m+ 2s/4>- m as 

I^WV-m = (-l) s det A«/det A^ . 
To calculate this explicitly, we use the following lemma. 

Lemma 4.1. For arbitrary q±, . . . , q n consider the matrix A = (ciy)j lJ '=i...n whose 
entries are = (qj) 1 — {qj)" 1 ■ Then 

det^ = n u ,2 i)' 2 - {^'\t - <k'\? 12 ) n (*■ - 



[<ii - ii 



i<j 



Expanding the determinant, it is easy to see that this formula is equivalent to 
the Weyl denominator formula for the C„-root system. There is also a simple direct 
way of proving it, using that q q Zq- 1 is polynomial in q + q" 1 . 

Applying the lemma, we calculate the determinants and find that 

lb ~ I -l^ tt ( 1/2 1/2 -1/2 -l/2\ / -1/2 1/2 1/2 -l/2\ 

w - m {q s -q s \Q> q/ -q« Qj ) [q° ij -is' q 3 1 ) 
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Substituting qj = q 



-m J r2j J rX 



we arrive after simple transformations at the formula 



m-\-2s 



n 



-m+J-l _ „m-J+l^ ^-rn+j-l+rr _ „m-j+l-a: 



3 



if>-m L \ (q ] - q~ j ) (q J+x - q- j ~ x ) 

In particular, for s = m we have : 

ipm /ip- 



n 

3=1 



q3 . q Jt^ 



(4.4) 



(4.5) 



Let us fix ip m in the following form: 



3 = 1 



Then the relation (|4.5|) gives us that 



n (^ + * - ■/ 



(4.6) 



3 = 1 



Note that our choice of ip m implies that all ip v will be Laurent polynomials in q x 
(see formula (4.4)). 

This allows us to prove the following proposition. 



Proposition 4.2. (i) The function ip(x,z) given by (4.1), (4.4) and (4.6) is a 
Baker- Akhiezer function for R — A\; 

(ii) tp satisfies the difference equation Lip — (q z + q~ z )ip , where L coincides with 
the operator (2.22): 



n x — m n —x+m 

L=- — ^— T x 



qX-\-m q—x—m 



T -i 



q^ — q * 

(Hi) ip is symmetric in x and z: ip(x, z) — ip{z, x). 



Proof. Part (i ) is p roven above. To prove (ii) we apply the standard argument due 
to Krichever | Krl |. Namely, let us consider the functi on <p = Lip — (q z + q~ z )ip . 
The first remark is that <f> still satisfies the conditions (4.2): (T| — T~ 3 ) (f> = for 
each j = 1, . . . , m and q 2z — 1. Indeed, it is obvious for Lip since the operator L 
does not involve z. Further, f(z) = q z + q~ z satisfies the conditions ( |4.2| ), hence, 
f(z)ip will satisfy them, too. 

Our second remark is that <p can be presented as follows: 



<P = cf 



3=0 



(p 3 q 



(-m-l+2j)z 



This follows directly from the formula ([LlJ) and our definition of (p. 

Let us calculate now the coefficient (po using ([0])and the definition of 
gives: 



(4.7) 



This 



q 



x-\-m 



-1p- m (x - 1) - 1p- m (x) = . 



q -q 

In the same way, (p m +i — 0. Thus, the expansion ( [4.7D contains m terms only, 
hence, it must be zero due to Lemma 3.2. This proves part (ii). 
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To prove part (iii), first notice that according to (4.4) and ( [4. 6| ) ip has no singu- 
larities in the cc-variable and it may be presented as 

m 

^ = 1 XZ E *»(*)<T' (4.8) 



We know that, a s a function of X, ip is an eigenfunction of the operator ( 2.22 ). 
Invoking Lemma |2.6| , we conclude that ip satisfies the following conditions in x: 

(Tl - T~ j ) ip = for each j = l,...,m and q 2x = 1 . (4.9) 

Thus, ip(x,z) must coincide with ip(z,x) up to a z-depending factor: tp(x,z) — 
F(z)ip(z 1 x). Switching x and z, we conclude that ip[z,x) = G(x)ip(x, z). This 
implies that F — G" 1 is constant. Expanding the coefficients ip v (x) with the help 



of the formulas (i.4) and (|Lq ), we see that ip(x,z) contains the term 

^ ^\m^-m(m+l)^cz^ma:+m2 

Since this term is symmetric in x and z, we conclude that F = 1 and ip(x, z) — 
ip{z,x). 

□ 

4.2. Normalized BA-function. Now we are going to extend the results of the 
previous subsection to the higher rank case. Above we have proved that a Baker- 
Akhiezer function ip is determined uniquely (up to a x-depending factor) by its 
properties ( |3.3[ )-( |3^5| ). Let us impose the following normalization condition on ip, 
prescribing its leading coefficient i/v> v to be the following: 

V = II lit? - («» *)] > [«] == <?" - • (4-10) 

aG-R+ j=l 

Definition . A n orm alized BA functio n is a (unique) function ip(x,z) with the 
properties (3.3)— (3.5) and normalization ( 4.1d| ). 



Let us discuss briefly the geometry of the polytope Af — supp^) defined by (3.4). 
We mentioned already that Af is the Minkowski sum of the segments associated with 
the vectors m a a v (a v € R+)- It is convenient to use a more symmetric definition 



of Af, which is obviously equivalent to (3.4) 



Af = {v = I ~\m a <L< \m a } . (4.11) 

aeR+ 

To better understand its structure, let us choose a generic direction v in V and 
consider the height function (v, •) on V. Then the highest and the lowest among 
the points of Af will be the points ^X)aei? + ^ n - m a<^ J , where the signs in the sum 
either coincide with the signs of (a v ,w), or are exactly the opposite. This shows 
that the vertices of Af have the form ^ . , s >Q r7i Q o; v , where v is a generic 
vector. For instance, taking v from the positive Weyl chamber C, i.e. such that 
(v, a) > for all a £ R+, we obtain p v as one of the vertices of Af. Other vertices 
will correspond to other chambers and will be of the form wp v with w € W. If 
we take now the vertex p v , then its adjacent vertices will correspond to the Weyl 
chambers, adjacent to the positive one. There are exactly n of them, n — ranki?, 
and they are the images SiC of C under the simple reflections s\,...,s n . It is known 
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that a simple reflection s, leaves the set i?+\a^ invariant, while sending on to —on. 
Thus, s.;/9 v = p v — rdiOtl ^ rrii := m Qi . By the way, this implies that (aj, p v ) = m„ 
i.e. p v can be rewritten in terms of the fundamental coweights bi as 

n 

/ 9 V =^m i 6 i , mi=m ai . (4-12) 
i=i 

So, the edges of TV coming out from the vertex p v are given by the vectors — m^a^ 
where a\, . . . ,a n are the simple roots from R + . In the same way, the vectors 
—rriiwa^ will lead from the vertex wp v to its adjacent vertices. Summarizing, we 
arrive at the following result. 



Lemma 4.3. The vertices of the polytope (4.11) have the form wp v , w € W. 
If v, v' are two adjacent vertices, then v 1 equals v + m a a v for a proper a G R. 
Moreover, in this case we have v' = s a v and (a, v') = —(a, v) — m a . 



Our next proposition justifies our choice of normalization (4.10) 



Proposition 4.4. The normalized BA function tp has the following properties: 

(i) for all w € W the coefficient ip u in the vertex v = wp v of the polytope M has 
the form 

^ P v= n lib' -fa*)]; ( 4 - 13 ) 

(ii) for all v lying on the (one- dimensional) edges of Af the corresponding i\) v are 
polynomial in x; 

(Hi) ip is quasipolynomial in both x and z. 

Proof. Let v and v' be two adjacent vertices of Af, so for a proper a £ R we have 
v 1 = s a v(i/) = v + m a a y . Let us introduce tpj — tl) u +j a v and qj = q 2 j- m °>+( a > x ) 
(j = 0, 1, . . . , m a ). Introduce also a matrix A with entries dj j = (qj) 1 — (qj)~ l (1 < 
i < m Q ,0 < j < m a ). In these notations, the conditions (3.5) can be expressed 



as Av — 0, where v is the column v — (tpo, ■ ■ ■ , ^ )r!a ). This system is completely 
analogous to the system ( |L3| ). Repeating the same arguments, we arrive at the 
formula 

^o^f[ [ - ma+J -fy ] ~^ +J - 1] , (4.14) 

where [a] denotes q a — q~ a . In particular, for s — m a we have the following 
expression for the ratio of the coefficients ip^ and i/v at two adjacent vertices of 
the polytope N: 

Recall now that v,v ! can be presented as wp v , w' p v for proper w, w' € W. 
Moreover, we have w' = s a w, and ±a are simple roots in the sets w'R + and wR+, 
respectively. This means that if we denote the common part of the sets wR+, w' R+ 



by S, then wR + = S U {-a} and w'R + = S U {a}. Now the formula fl4.13| ) in 
question gives us that 

/ / / rf b' - ("> x )\ 



3- 
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A [j + (a, x)] 



Since this expression coincides with the formula (4.15), and since ( 4.13Q is valid at 



one vertex v = p v , we conclude that it must be valid at all other vertices of Af. 
This proves part (i) of the proposition. 



The second part follows directly from the formulas (4.13) and (4.14). To prove 



part (iii), we have to show that all coefficients ip v are polynomial in x. The BA 



function in formula (3.16) is clearly quasipolynomial in x, so we only have to prove 
the absence of singularities in the normalized BA function ip. 

Suppose ip has a pole along a hypersurface F(x) — 0. Without loss of generality, 
we may assume that it is irreducible. Multiplying ip by F{x) we will obtain the 



function if) which will be quasipolynomial in x sharing with if) the properties (3.3) 
(fT5|). Now take generic point xq on this hypersurface, so F(xq) = 0. Then all the 
coefficients ip v = F(x)ip v for v lying on the edges of the polytope Af must vanish at 
x — xo (because ip u are polynomial, see part (h)). However, ip does not vanish at 
Xq (otherwise if) would be nonsingular at F = 0). Let us denote by Af the support 
of if)(xo, z) (i.e. the convex hull of all v with nonzero ip v (xo)). We have, therefore, 
a strict inclusion 

Af^Af. (4.16) 

Moreover, this shows that the polytope Af has no vertices lying on the edges of Af. 
This will lead us to a contradiction in a moment. 

Indeed, take any a € R and consider the function f(x) = q 2 ( a > x > restricted to 
the hypersurface H : F = 0. Since H is irreducible, / either is constant along H or 
takes infinitely many different values. Suppose now that for every a £ R the latter 
alternative hold s. T hen we would be able to choose generic xq € H such that our 



proof of lemma 3.2 would work, i.e. such that all m a x m a determinants arising 



from considering q-s trin gs inside the polytope Af would be nonz ero. H owever, this 



would imply lemma 3^ and this would contradict the inclusion ( 4.16 ) . 

The upshot is that for some a € R the function g 2 ^- 21 ) is constant along our 
hypersurface H. But in this case we can choose generic xo £ H such that lemma 
is still applicable to all roots except a. Hence, the support Af of i/> due to lemma 
must contain a polytope Af' , which by definition is the Minkowski sum of all 



3.2 



3.3 



m/3/3 v («^/3g R+)- However, Af = Af'#m a a v . This means that each copy of Af' 
lying inside Af, has at least one vertex on the edges of Af, but Af has no vertices on 
the edges of Af. This contradiction proves that if) has no singularities, therefore, it 
is quasipolynomial in x. □ 



So, our choice of normalization (4.10) guarantees that ip will be quasipolynomial 



in x. In a certain sense, it is the "minimal" quasipolynomial BA function. 

Corollary 4.5. Any BA function ifj(x, z) which is quasipolynomial in x is the 
normalized BA function multiplied by some polynomial in x. 



This follows immediately from the formula (4.13) which shows that all the coef- 
ficients ip wp v have no common divisors (as polynomials in x). 

4.3. Duality. A remarkable and important property of the normalized Baker- 
Akhiezer function is a certain duality between x and z variables. Namely, it turns 
out that switching x and z leads to the BA function for the dual root system. In 
particular, for A,D,E root systems we have simply ip(x,z) = if>(z,x). The next 
proposition is the main step in establishing this symmetry. 
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Proposition 4.6. The normalized BA junction ip satisfies the following conditions 
in x: for each a G i?+ and s = 1, . . . , m a 

ip(x + ^sa v ,z^J = ip(x-~sa v ,Z^ for q 2 ^ = 1 . (4.17) 



Proof. Take any a € R and consider two adjacent vertices v, v' of the poly tope M 
such that v' = v + m a a y . As we know, v and v' can be presented as u>p v , w' p y 
for proper ui, u/ € VK. Moreover, denoting the set w' R+ U u>i?+ by S, we will have: 

= SU{-a}, w'R + = S(J{a} 1 s a (S) = S . 



Consider the function ip — I T x 2 
quasipolynomial in z, 

i> = q 



ip. Similar to ip, it will be 



2(x,z) 



tEP x 



a T q 



2(t,z) 



(4.18) 



and it is clear that its support lies inside the polytope 

Kf = Af'#l , 



isa v and 



(4.19) 
v v and W', in 



where I is the segment with the endpoints at v 

its turn, is the Minkowski sum of all mp(3^ with (3 € S (we recommend the reader 
to draw a picture for A2 case). 

Now take xq such that g 2 ( a > x °) = 1. We would like to prove that for such a xo 
ip will be zero. Suppose that it is not the case. We may assume that Xq is generic 
enough, i.e. such that we still can apply Lemma 3.2 to all (3 E S. This, together 
with Lemma 3.3, implies that supp^) can be presented as Af'fptt for some domain 



f2. Comparing with (4.19) gives us that f2 C I. This implies that at least one of 
the coefficients a T in (4.181) with r € I must be nonzero. So, we will arrive at a 



contradiction as soon as we prove that all a T for r G I vanish. We will see in a 
moment that this essentially reduces the problem to the rank-one case. 

Indeed, notice that q 2 ^'^ ^2 Te ia T q 2 ( T ' z \ by its construction, coincides with 



ip( ', where 
^ (0) := 



2(:r, z) 



2<r,z) 



is the part of ip, corresponding to the edge \v, v'\ of the polytope N '. We recall that 
we consider x — x such that q 2 ^ a - x ^> = 1. 
So, we need to show that 



')^ (0) =0 for q 2 ^ = 1 . 



(4.20) 



If our ip^ were the normalized BA functi on for the root system R^°> = {±a}, 
this woul d follo w dire ctly from Proposition |4.2| (ii). This is almost the the 
formulas ( 4.14| )-(4.15) clearly indicate. The only difference is that, due to (4.13), 
the coefficient \p v with v — wp y looks as follows: 



^=]ll[\j-(fi,x)]l[[{a,x)+j] 

PeSj=i j=i 
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(here we used that wR+ = S U {—a}). So, the difference with the similar formula 



(1.6) for the rank one case comes from the factor rises Ilfcib — (0i x )]- However 



this factor is invariant under reflection s a , hence, it does not affect the properties 



(4.2C). □ 



Next thing to prove is that the support of ip in the ir-variable is the following 
polytope A/" v : 

7V V = {v = p - l a a\0<l a < m a } . (4.21) 

a£R + 

To prove this we look first at the formula ( |3.16 ), which gives an expression for 



another BA function with ip p v normalized as in (3.17). It is clear from (3.16) that 
such ip is also quasipolynomial in x with supp^ = supp x ip p v . If we renormalizc 
now ip in order to get ip p v as in (4.1C), then we will still have that the support of ip 



equals the support of ip p v which now will coincide with 7V V . So, for the normalized 
BA function ip(x,z) one has snpp x ip = Af w . 

As a result, we see that the normalized BA function can be presented in a form 

with the summation taken over \i £ J\f v and v £ J\f, respectively. The highest term 
in this expression will correspond to p = p, v = p v . 

Summarizing, we see that the properties of ip(x,z) in the x-variable are com- 
pletely analogous to its properties in z. Thus, we obtain the following important 
property of the normalized B A function which reflects the symmetry between x and 
z variables. 

Theorem 4.7 (Duality). Let ip(x, z) denote the normalized Baker- Akhiezer func- 
tion associated to a root system R with multiplicities m — {m a }, and ip v (x,z) 
denote a similar function associated to the dual system R v with m a v = m a . Then 
ip{x, z) = 4^ '( z i %)■ 



Corollary 4.8. Let ip be the normalized BA function for a given root system R and 
m = {m a }. Let to € P and ir G P v be (quasi) minuscule weight and coweight for 
the root system R. Then ip(x, z) solves the following bispectral system of difference 
equations: 



Here D x and are the corresponding Macdonald operators ( |2.2| )-(2.3) and ( |3.7|) 



( p.q ), while m\ stands for the orbitsum (2.8) 



5. Algebraic integrability and applications 

5.1. Algebraic integrability. Above we have shown that the normalized Baker- 
Akhiezer function ip(x, z) associated to a datum {R, m} is an eigenfunction of the 
Macdonald operators D x , where it is any (quasi)minuscule coweight for the root 
system R. In fact, ip is a common eigenfunction of a much bigger commutative ring 
of difference operators. This follows in a standard way from its analytic properties in 



the z-variable (cf. [Krl|). To formulate the result, let us recall the definition of the 



ring 5H V which consists of all polynomials f(z) of the form f(z) = X^eP v fv<f^ v ' z . 
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satisfying conditions (B.£) for each a £ R and s = 1, . . . , m Q . It is easy to see that all 
W- invariant polynomials belong to the ring 9\ v . According to Chevalley theorem, 
W^-invariants form a polynomial ring with n = ranki? generators , . . . , m& n , 
which are the orbitsums for the fundamental coweights b±, . . . ,b n . However, the 
ring 9l v is much bigger, for instance, it contains the principal ideal generated by 



the polynomial Q given by (3.13) 



Theorem 5.1 (Algebraic integrability). For each polynomial f(z) from the ring 
5H V there exists a difference operator Df in x on the lattice P v such that Dfty = 
f(z)^!. All these operators commute. For a (quasi) minuscule coweight n of R and 
f = m 1T (z) t he c or resp onding operator Df coincides with the Macdonald operator 



given in d2.l|)-(p.3j) 



Proof. Everything is based on the following result. 

Lemma 5.2. Any quasipolynomial in z of the form f> — q 2 ( x - z *> ^2 vepv <j) u (x)q 2 ^- z ^ 
which satisfies the conditions ( |3.5| ) can be obtained by applying a proper difference 
operator in x to the BA function ipix, z). 



To prove the lemma, we recall that according to lemmas B.2 and 3.3 the support 
of <f) can be presented as No#Af, where Af is the support of i/j and Nq is some convex 
polytope. Now choose any vertex v of the supp</> and let r be the corresponding 
vertex of Nq such that v £ r + M C supp0. We can "kill" the coefficient (\> v of <f> 
by subtracting a function Tf.if taken with a proper coefficient d{x). The resulting 
function (f> = (j) — d(x)Tf,ip still satisfies the conditions of the lemma, but it has a 
smaller support. Repeating this, we will eventually get zero, and this proves the 
lemma. 

Now, to prove the theorem, we notice that for / € 9^ v the function <fi = f(z)ip 
will satisfy the conditions of the lemma, hence, <f> = Df?p for a proper difference 
operator Df — J2veP v dvT£. All these operators commute since they have ip as 
their eigenfunction. Indeed, for /, g £ £H V we have [Lf,L g ]ip = (fg — gf)i/j = 
0. However, one shows easily that if a difference operator M in i annihilates 
a (nonzero) quasipolynomial if) in z, then M = (otherwise suppM-0 would be 
nonempty). Hence, [Lf,L g ] = 0. 

Finally, for f(z) = m 7r (z) and the corresponding Df we will have Dfip = nv(z)i/>. 
The same is true for the Macdonald operator D^.: D^ip = m 7V (z)t/j. Hence, these 
two operators coincide. □ 



Remark. The duality between x and z implies that the normalized BA function will 
be also a common eigenfunction of a " dual" commutative ring of difference operators 



in z variable, isomorphic to the ring 9\ of polynomials with the properties (2.15) 
Thus, we have a bispectral pair of commutative rings, in the spirit of jW|]. 



Remark. In the limit q — > 1 this theorem reduces to the result from CSV | , where 



the algebraic integrability was established for the quantum trigonometric Calogero- 



Suthcrland-Moser problem and, more generally, for its generalizations |OP| related 



to the root systems. Proof in [CSV] was based on results by Heckman and Opdam 
|HO| , ||], who developed a nice theory of multivariable hypergeometric functions 
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related to root systems. An independent direct proof was obtained in [Chi]. For 



A„_i case, i.e. for the Macdonald-Ruijsenaars operators (2.7) with t = q m (to 6 Z) 
the algebraic integrability was established by Etingof-Styrkas in ]ES| . They con- 
structed ip in terms of certain intertwiners from representation theory for quantum 
groups. Their formula for ip was made more explicit by Felder and Varchenko in 

IfT 



5.2. Liouville integrability. As a part of the previous theorem, we have con- 
structed n = ranki? commuting difference operators D; corresponding to the ba- 



sic IF-invariants / = m^z), which are the orbitsums (2.8) for the fundamental 



coweights hi € P v . Those of 6j which are (quasi)minusculc, will lead to the Mac- 



donald operators (2.1), (2.2) with t a — q~ ma . For others, there is no such simple 
explicit formula. However, it is easy to evaluate leading terms in Di and they look 
similar to the leading terms in Macdonald operators. 

Proposition 5.3. For any coweight tt G P v and its orbitsum f = m 7r (z) the 
operator D j constructed in Theorem \5. J| will have the form 

Df = £ ".- /" + • • • , 

t£Wtt 

where the dots stand for a sum of the terms a v T v with v € tt + Q v lying inside the 



convex hull of the orbit Wir, and the leading coefficients a T are given by (2.5), ( |2.6| 
with k a = —m a . The operator Df is W -invariant. 



Proof. Formula for a T follows directly from the construction of Df in theorem ( p.l 



and formula ( 4.13 ) for the leading coefficients in if). It remains to prove that Df is 



IF-invariant. To this end we have the following symmetry property of if) . 

Lemma 5.4. The normalized BA function is VK-invariant in the following sense: 
tp(wx,wz) — tp(x,z) for any w £ W. 



To prove this we notice that for the leading coefficients (4.13) of ip one nas 



i pwp v ( wx) — ip p v(x). Since ip(x,z) :— ip(wx,wz) shares with ip the properties 



( p.3| )— (3.5), they must coincide due to the uniqueness of ip. 

Using the lemma and VF-invariance of /(z), one gets that [Dip](wx,wz) = 
f(wz)ip(wx 1 wz) = [Di/)](x, z). From this it easily follows the property a wv {wx) = 
a v {x) for the coefficients of D, i.e. its W^-invariance. □ 



Let us look what happens if we change the multiplicities m a . The formula ( p.6[ ) 
shows that the leading coefficient in Df will be rational in t = {t a }. In fact, it 
is not difficult to prove that all the coefficients of the operator Df with / = rrv(;z) 
will be rational in t. 

Lemma 5.5. All the operators Df depend rationally on t. 



Proof. Recall that the construction of Df was given in Lemma 5.2. Now let f{z) — 
m 7! {z) be an orbitsum, tt 6 P + . We can construct Df as in Lemma |5~2| , starting from 
<f) = f(z)ip and "killing" at each step the highest coefficient of cf>. This shows that to 
calculate Df in this case it is sufficient to know a fixed number of the coefficients tp v 
in the normalized BA function ip, namely, those with v — p v € — tt + Wtt. Thus, if 
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all these coefficients were rational in t, then Df would be rational in t, too. This is 
not the case, however, because already ip p v is not rational in t. Let us renormalize 
ip in such a way that ip p v = 1, and let us calculate Df using this ip. Of course, the 
resulting operator Df will differ from Df, but the relation is simple: 

D f = AoDfO A" 1 , 

where A is 

A= II H(q- {a ' x} -q- 2l+{a ' x) Y , (5.1) 
a£R + j=l 

which is a specialization of ( |2.5| ) in case t = q~ m . An important point is that 
despite the fact that A is not rational in t, a ratio 

A(x)/A(x + v) 

is rational in t for any v 6 P v . Therefore, it is sufficient to prove that all the 
coefficients of the operator Df depend rationally on t. 

Summarizing, we showed that it is enough to check that for any fixed r € 

the coefficient ip v with v = p v — t will be rational in t. Recall now that the 
renormalized ip, as well as ip itself, satisfies a difference equation in the z-variable: 



derived in Theorem 3.7. Using this equation, one can calculate the coefficients of ip 



recursively, starting from ip p v (one gets this recursion similar to [HO], expanding 
the coefficients of the operator in a series in z). Since the Macdonald operator 
is polynomial in t, each ip v with v = p v — r, calculated from this recursion, will be 
polynomial in t. This completes the proof. □ 

As a corollary, we obtain that all the operators Di related to the fundamental 
coweights, must commute for all values of t (since they commute for t = q~ m with 
any integer to). Let us consider one of the Macdonald operators as a Hamiltonian 
of the corresponding quantum problem, which can be viewed as a generalization 
of the trigonometric Ruijsenaars problem to all root systems. Then what we just 
proved is the complete (Liouville) integrability of this quantum problem. So, in this 



way we recover the result obtained by Cherednik [CI] 



Corollary 5.6. A quantum problem related to an arbitrary Macdonald operator is 
completely integrable. 



5.3. Shift operators. Another important result by Cherednik is the construction 
of the so-called shift operators, which are g-versions of the operators constructed 
by Opdam for the case q — 1. Let us explain how they appear in our approach. 

Let us consider a root system R and two sets of integer multiplicities m, m. 
Suppose that m > m, i.e. m a > fh a for all a € R. Denote by ip,ip corresponding 
normalized BA functions. Comparing their properties (3.5), we then apply Lemma 
5.2 to conclude that there exists a difference operator S + in x such that ip = S + ip. 
Let D, D be the Macdonald operators (2.2) with t = q~ m , i = q~ m , so D^tp = m^ip 
and D w -ip = m^tp. Hence, 

(DoS + -S + o D)4> = Dip- S + {m^) = . 
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Therefore, we arrive at the identity 



DoS + = S + oD. 

It means that the operator S + intertwines two Macdonald operators with different 
t. In the same way, S + intertwines each pair Di , Di of the corresponding operators 



from Proposition 5.3 



Consider now the case t = qt (i.e. fh a — m a — 1 for all a). We will use subscripts 
denoting by D m the Macdonald operator with t = q~ m . So, we have shown that 
for any integer m = {m a } we have two Macdonald operators D m , D m -\ and their 
intertwiner with the relation: 

D m o S+ = S+ o • (5.2) 



Similar to Lemma 5.5, one proves that the intertwiner 5+ will depend algebraically 
on t = q~ m ', so the relation (5.2) will make sense for any value of to. Corresponding 
is called a shift operator, since it sends eigenf unctions of D m _i to eigenfunctions 
of D m . From its construction we can calculate easily the leading terms in S m . 
Introduce g v G P v as 



\ E a " ■ < 5 ' 3 ' 



Proposition 5.7. The shift operator above has the form 

S+= E b T T T + ..., 

where the dots stand for a sum of lower terms b v T v with v G g v + Q V lying inside the 
convex hull of the orbit Wg w . The leading coefficients b T are given by the formulas 

bgV {x) = q E a>0 m a Ai- m (^ + e v ) b {x) = b {w -i x) 
A_ TO (x) 



where g is given by (5^) and A_ m stands for the function (5.1). The operator 
is W -invariant and it sends the normalized BA function ip m —i to ip m . 

For a proof one should compare the supports and the leading coefficients in two 
BA functions, tp m and ip m -i. 

Another shift operator will appear if we consider, as above, two normalized BA 
functions ipm,?Pm-i corresponding to to = {to q } and fh = m — 1, and multiply 
iptn—l by a polynomial c m (z)c m (— z) where c m (z) is defined by the formula: 

c m {z) = [] K " (« v > z )l - H : = i a - q~ a ■ (5.4) 

aeR+ 

It is easy to see that the resulting function <p = c m {z)c m {—z)ijj m ^i will satisfy the 
conditions (3.5). Hence, due to Lemma |5.2| , it has the form S~ip m for a proper 
difference operator S~ = S"~ (in the a;- variable). Its leading terms, again, can be 
calculated easily. Similarly, we get the intertwining relation D m _i o S 1 " = S 1 " o D m 
which, again, extends analytically to all values of to. Next proposition summarizes 
the properties of the constructed S~ . 
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Proposition 5.8. For integer m = {m a } there exists a difference operator S — 
S m which intertwines the Macdonald operators D rn , £) m _i: 

D m -l o S rn = S m O D m . 

It has the form S~ = J2 T eWg v b T T T + . . . , where the dots stand for a sum of lower 
terms b v T v with v G g v + Q v lying inside the convex hull of the orbit Wg v . The 
leading coefficients b T are given by the formulas 



b e ,(x) = (-i)l^l g -S«>.". ^fr + g-} , b we ,(x)=b e ,(w- 1 x) 

Ai- m {x) 



in accordance with (5.3), ( p.l| ) (here \R+\ denotes the total number of positive 
roots). Application of S m to the normalized BA function ip m gives S m tp m = 
c m (z)c m (—z)ip m _i, where c m is defined in (5.4). 



Remark. In case when R consists of two orbits, R = R\ U R2, one can apply 
similar arguments for the case when multiplicities increase for one of the orbits, i.e. 
m = rh + lj with 1, being the characteristic function of Ri C R. This proves the 
existence of the intcrtwiners, shifting from mtomllj. Corresponding versions of 
propositions 5.7, 5.8 are straightforward. 



Remark. More generally, we can consider a pair of data {R, rh} and {R, m} with 
R C R and with fh a < m a for all a £ R. This leads to the intertwiners between 
Macdonald operators for different root systems. Unlike the s hift operators S m 
above, they exist for integer values of m only. As an example, see | CV where such an 
intertwiner was constructed explicitly for the case q = 1, and R = A n C R = A n+ i. 



We conclude the section by a description of all possible shift operators. First we 



need a result showing that the commutative ring from Theorem 5.1 is maximal in 
a certain sense. 

Proposition 5.9. Let D = be a Macdonald operator related to a root system 
R and t — q~ m with m a G Z + and let 5H V be the ring of polynomials in z with the 
properties ( |3.9[ ). Suppose that a difference operator L commutes with D and has 
rational coefficients. Then L is one of the operators constructed in Theorem 5.1, 
i.e. L — Df for some f G 9\ v . 

Proof. Suppose that [L, D] = and consider a function $ = Lip where ip is the 
normalized BA function. Under assumptions of the theorem, we have that <fr(x, z) 
is of the form $ = q 2 ( x ' zS >tp where <p is polynomial in z and a rational function in 
x. From the commutativity of D and L we get that $ is also an eigenfunction for 
the Macdonald operator D: £>$ = m % (z)^. It is not difficult to deduce that in 
case t a = g~ m ° with m a G Z+ any eigenfunction of D = either has infinitely 
many poles along hyperplanes of the form q^ a ^ = const, or has no poles at all. 
By our assumptions, $ cannot have an infinite number of poles, hence, it must be 
quasipolynomial in x. Then from the same equation, using lemma B.6, we derive 
that suppj.^ must coincide with the polytope (3.4). Further, with the help of 
lemmas 2.6, 2.7 we conclude that $ satisfies the conditions (4.17). Hence, $(z,x) 
is a Baker- Akhiezer function for the dual root system R v . Now corollary [4. 5| implies 
that $(z,x) — f (x)i/j v (x , z) , where ip v (x,z) is the normalized BA function for i? v 
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and / is some polynomial. Switching x, z and using duality, we obtain that $ must 
be proportional to the normalized BA function: $(x, z) = f(z)ij)(x, z). Notice that 
since 4> was obtained from ip applying a difference operator in n, it shares with ip 



the properties (3.5) in the z-variable. Since this must be valid for any x, we obtain 



that / itself must be from the ring (3.9). Thus, / belongs to 9\ v and L = Df. □ 



Proposition 5.10. For positive integer to — {m a } let D m and D m _i be two 
Macdonald operators with t = q~ m and t = q~ m+ , respectively. Let S be any 
difference operator in x with rational coefficients satisfying the intertwining relation 
D m oS = SoD m —i. Then S = S + oL where S + is the shift operator from Proposition 



5.1 and L commutes with D m -\. 

Proof. Let ip m and ip m -i be the normalized BA functions for multiplicities to and 
m = m— 1, respectively. From the intertwining relation we obtain that ip := Stp m -i 
will satisfy the same difference equation as does ip m : D m ip = m 7r (z)jp. Due to 



Proposition 5.E , ip must have the form f(z)ip m for some polynomial f(z). The same 
arguments as before give that f(z) must belong to the ring £K V related to to = m—1. 
Hence, Sij) m -i = (S + o L)ip m —i where L = Df is such that L"0 m _i = f(z)ip m -i. 
This implies that S = S + o L. 

□ 

5.4. Relation to Macdonald polynomials. Let ip(x, z) be the normalized BA 
function constructed as above starting from the data {R, to}. Let us consider two 
functions <i>± obtained from ip by (anti)symmetrization in x: 

= ^ wx > z ) ' *- = H {-l) w i/j{wx, z) . 
Notice that the same will be the result of (anti)symmetrization in z, due to Lemma 



5.4. So, the resulting functions <£>± are (anti) symmetric in z, too. In fact, the func- 



tions <I>± (x, z) are closely related (for special z £ P + ) to the Macdonald polynomials 



P\(x;q,t) . Let 8 — S rn (x) be the function defined by ( 2.30 ), that is, 



s(x)= n n 

qG-R-j. j=—mot 



Recall also the definition (3.1) of p = p m : 

p m = ^ ^ m«« ■ (5-5) 

In particular, in case to = 1 the vector p\ coincides with g from (p^). 

Theorem 5.11. Let us substitute z = A £ P + into (f>±. Then 
$+ {x, A) = c(A, m)P x+Pm {x; q, q~ rn ) 

and 

(x, A) = c(A, m)8{x)P x . Pm+1 (x; q, q m+1 ) , 
where 5 is defined above and the factor c(A, m) is given by the formula 

c(A,m)= [] I[b'-( aV ' A )]- ( 5 - 6 ) 

aeR+j=l 
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We suppose that (A, a v ) > m a for all a G R + to ensure that c(A, m) 7^ 0. In other 
words, A G /O m +i + P+. 

Proof. Let 7r be any (quasi)minuscule coweight of R and D — denote the 
corresponding Macdonald operator ( |2.2| ) with t = g~ m . Each of the functions 
ip>(wx, A) = ip(x, w~ 1 \) satisfies the same equation Dip = m„.(A)'0, hence we have: 

D$±=tiw(A)#± m*(A) = 9 2(T ' A) ■ (5-7) 



Using Lemma 2.9, we conclude that = <I> + and p(x) = S 1< £'_ will satisfy the 
equations 

Dp^m 7r (X)p, L»p = m 7r (A)p, 



where D and D are the Macdonald operators (2.2) with t = g _m and t — q m+1 , 
respectively. These are exactly the defining equations ( 2.13 ) for the Macdonald 
polynomials P\+ Pm (x; q, q~ m ) and P\- Pm+1 (x; q, q m+1 ). So, to prove the theorem, 
we only have to check that p and p are symmetric polynomials and to calculate 
their leading terms. 

The VF-invariance of p,p is obvious since $± is (anti)invariant. Calculating the 
leading term is quite straightforward since we know the leading terms in ip(x,\). 
So the only non-trivial thing to prove is that is divisible by S(x). Here we can 



use the properties (4.17) of the BA function. Let us rewrite them after shifting in 
x by — isa v : 

ip (x, z) - if> (x - ja v , z) — for q 2 ^'^ = q 2j , j = 1, . . . , m a . 

Choose now any a £ R and x such that (a,x) = j, then its image x' — s a x under 
reflection s a will be x' — x~ ja v , so according to the property above, we will have: 

ip(x, z) — if){s a x, z) = if{x, z) — ip(x, s a z) = for (a, x) = j . 

Now let us split the sum $>-(x, z) = ^2 weW {—^) w if{x, wz) into pairs of terms with 
w, w' — s a w. As a result, we obtain that 

$_ (x, z) = for (a, x) = 1, . . . m a . 

<f>_ also vanishes for (a, x) = due to its antiinvariance. Invoking all a S R, we 
arrive at the following result. 

Lemma 5.12. For any a G R+ the function $_(x,z) = Xitogn/( _ ^) l "^( 2; ' wz ) 
vanishes along the hyperplanes (a, x) = 0, ±1, . . . , ±m Q . 

This is true for any z. However, for z = A G P the functions <i>± are quasiperiodic 
with respect to the lattice C = ojP v where u> = iri^ogq)^ 1 . Indeed, under a shift 
x — > x + / (I G C) each of ip(x,wX) (w G W) gets the same factor as does the 
function q^+P^f These translation properties imply that as soon as <E>_ vanishes 
for 2(a,x) = 2j, it will also vanish for q 2 ( a ' x '> = q 2 i . This proves that $_ for 
z = A G P is divisible by the polynomial 5 (x), thus, completing the proof of the 
theorem. 

□ 



Remark. The proof shows that $-(x, A) will be zero for all A G P+ which do not 
belong to p m +i + P+. Indeed, in this case must be divisible by 6(x) but has a 
smaller support. 
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Our expression for P\ +Pm (x;q,q~ m ) via the function <f> + is well-defined for suf- 
ficiently large A, namely for A € p m +i + P+- For smaller A it doesn't work, which 
reflects the known fact that in case t = q~ m with m a S Z + some of P M are not well- 
defined. Let us call the Macdonald polynomials P p (x;q,q~ m ) with (//, a) > 2m a 
for all a > the stable Macdonald polynomials. They are always well-defined (if q 
is not a root of unity). The following "localization" property is a direct corollary 



of Theorem 5.11 



Corollary 5.13. For t = q m with m a G the stab le Macdonald polynomials 
P\(x;q,q~ m ) are localized, i.e. in their expression ( 2.1 0| ) through the orbitsums 



P\ = m x + a ^ x 



only the terms with v = A — X^ae-R ^ a0t an< ^ ^ a = 0, 1, . . . , m a can appear, all other 
coefficients a\ v vanish. As a result, the total number of nonzero terms remains 
bounded as A increases. 



Below is a simple example which illustrates Theorem 5. If 



Example. Let us consider the case R — A± with m = 1. Using the results of 



section 4.1, we obtain (switching x, z) the following formula for the normalized BA 
function: 

lf)(x, z) = (q 1+z - q- l - z ) q^ x + (g 1 -* - «T 1+Z ) q^ x . 
As a result, we get the following expressions for Q± = ip(x, A) ± f>(—x, A): 

*±(x, A) = (q 1+X - q- 1 -*) (q^ x ± q^ 1 ^ 

+ (<? 1_A - q- 1+X ) (q {X+1)x ± <Z ( ^ 1)a; ) ■ 
Taking integer A > 2, we obtain according to Theorem |5.1l| the following formulas: 



P x+1 (x;q,q- 1 ) = (q 1 - X -q- 1 + X ) 1 $+ = m A+1 + |^jm A -i 
P\- 2 {x;q,q ) - 



[I - \}[x - l}[x}[x + 1} 



Remark. In the case 



= the BA function is a pure exponential, tp 
Theorem 5.11 reduces in this case to Weyl's character formula. Thus, one can 
think of Theorem 5.1l| as a generalized Weyl formula with ip being a " perturbed" 
exponent. Note that there is a similar result involving "non-symmetric" Macdonald 
polynomials (see [M3Q. Our ip(x,A) for A e P is polynomial, too, but it differs 
from non-symmetric Macdonald polynomials. In fact, if) is related to a natural q- 
analogue of the hypergeometric function by Heckman and Opdam [HO|. However, 
we will not discuss this relation here. 



Remark. In case R = A n our formula for <!>_ is the formula conjectured by Felder 
and Varchenko [FV] and proved by Etingof and Styrkas [ES] in the context of the 
representation theory for quantum groups. 
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The following proposition is another direct corollary of Theorem 



5.11 



Proposition 5.14. Let P^™ 1 ' = P tl (x;q,q m ) denote the Macdonald polynomial for 
a root system R and t = q m , where m = {m a } are positive integers. Consider the 



shift operators SiE constructed in Propositions 5.7 and 5.t. Define S^z as 



$ m -i o S m o S„ 



Then these operators act in the following simple way onto Macdonald polynomials: 



q+ p( m ) 
+ p( m ) 



c m (A)Pj; 



(to+1) 



-p m +i 



for all A G Pm+i + P+ arid S^P^_ = otherwise. Further, 



A- p (m+l) _ (_\\pV m ) 
D m- r X- Pm+1 — c m{ A )^\- p 



(m) 



for all A £ Pm+i + P+- In these formulas c m denotes the function (5.4). 



Proof. Let ip m (x,z) denote the normalized BA function for the system R with 
multiplicities m. According to Theorem |5. 11 we have: 



(m) 



(x,X) = ]T (-l) w ip m (x,w\) = c(X,m)S m P { x m +y +i (x) 



w£W 



As we know, the shift operator S m sends ip m (x,wX) to c m (A)c m (— \)4> m -i{x, wX) 
(here we use W-invariance of c m (z)c m (— z)). Applying it to we obtain: 

S-*™(x,\) = c m (X)c m (-X)^ m) . 
After rewriting it in terms of Macdonald polynomials, we get: 



S ,, 



c(X,m)S m Pi n :^ ] = c m (X)c m (-X)c(X,m- l^-iP^, 



which leads directly to the last formula from the proposition. Another part, involv- 
ing S+, can obtained in the same way. □ 



Later we will apply these results to prove the so-called norm identity for Mac- 
donald polynomials. In order to do this, we need one more result. 



Proposition 5.15. Let (,)k denote the scalar product ( 2.12| ), depending on the 



parameters k a which are supposed to be positive integers. Consider two shift op- 
erators Si , constructed in the previous proposition. Then S^ is "adjoint" to S£ 
in the following sense: 

<s+/,. 9 ) fc+1 = (-i^+'^-N/A^k- 

for any two polynomials /, g. 
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Proof. First, recall that according to [Ml| each of the Macdonald operators D k , 
being restricted to the space of VT-invariant polynomials, is self-adjoint with respect 
to the scalar product (,)/.. We should warn the reader that this is not true for 
k a G Z_. Nevertheless, for all k the Macdonald operator D k is formally self-adjoint. 
This is a purely algebraic statement, related to the following rule of calculating the 
adjoint: 

(a T T T )* = (^k\) ^ oT T o (A k A k a T ) , (5.8) 



where / = f{—x). The point is that this definition leads to the relation 

(a T T T f,g) k = (f, {a T T T )*g) k 

as soon as we are sure that a T AfeA fe is polynomial. 

In a similar way, let us define another type of adjoint by the rule 

(b v T»y = (A fe A fe ) _1 ofo (A k+1 A k+1 b u ) . (5.9) 
This definition leads to the relation 

{b v T v f,g) k+1 = {f,{b v T v )*g) k 

as soon as we know that b v A k+ iA k+1 is polynomial. 

Now let us consider the shift operator S k from Proposition 5.14 . As we know, 
it satisfies the intertwining relation 

D k +iS k = S^D k , 

where D k now stands for th e M acdonald operator with t = q k . Taking the adjoints 
according to the definition (5.9), we obtain that 



D k +\ — D k (S k 



(here we used that D k and D k+ i are formally self-adjoint in the sense of (5.8)). 

Notice that this coincides with the intertwining relation for the shift operator 
ST. One easily compares their leading terms, using the formulas from Propositions 
5.7 and |5.8| and the definition (5.£) of the adjoint (in doing this it is useful to 
present the function ( [2.30 ) in terms of ( pT5| ) as S m = (— 1)1-^+1 A m+1 /A_ m ). A 
straightforward calculation gives that the leading terms coincide up to a factor 
(— l)l- R +lg 2 S k " . Thus, the intertwiners should also coincide, due to Proposition 
5 . 10| . This proves that 

(S£r = (-i)i*vE*-§r, 

where the adjoint is understood in the formal sense (5.9). 

Now, in order to derive the relation between the scalar products, it remains to 
check that if we present the operator S k as a sum ^2 o v T v then each coefficient b v 
after multiplying by Afe + iA fe+1 becomes polynomial. Clearly, it is enough to prove 
that the operator A — o <5 m _i has polynomial coefficients. To prove this, we 
consider a quasipolynomial (p(x,z) — q 2 ( x ' z ' 8 m -\(x) which obviously satisfies the 
conditions (4.17) for j = 1, . . . , m a — 1. By Lemma |5.2| , can be presented in the 
form (f> = L z (ip m -i) for some difference operator L in the z-variable, where VVn-i is 
the normalized BA function. Using that S^) m -\ = ip m , we see that S^fj) = L z ip m , 
i.e. it is polynomial in x. This proves that the operator A = 5+ o 8 m —x maps 
polynomials into polynomials, thus, it must have polynomial coefficients. This 
completes the proof of the proposition. □ 
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5.5. Macdonald Cherednik identities. Our results imply immediately three 
Macdonald's conjectures about polynomials P\, which were first proved (for all 
reduced root systems) by Cherednik JCl| , ^2| |. 

First, let us derive the norm identity in a traditional way, following the idea 
first used by Opdam Jo| in case q = 1 . Let , as above, denote the Macdonald 
polynomial Pfi(x; q, q k ) related to the system R and parameters k a S Z+, and (, )k 
stands for the corresponding scalar product (2.12). Now we take an arbitrary weight 
A G Pm+i + P+ and rewrite the scalar product (PxH^ +1 , Px^p^ +1 )m+i with the 
help of propositions 5.14 and 5.15| : 

/p(m+l) p (m+l) \ _r / 5. / c+ p( m ) c+ p( m ) \ _ 

(-l)l^l^-[ Cm (A)]- 2 (pg m ,5-5+P^ m ) m = 

This reduces the problem of calculating the norms at t = q m+1 to the same problem 
but for t — q m . For instance, in the simplest case when all m a are equal, we may 
descend to m — 0. Macdonald polynomials P x are simply the orbitsums mj(i), 
so in this case (P x °\ P^)q — \W\. This gives the norm identity: 

aeR+ 3 = 1 > A > J\ 

This result extends easily to the two-orbit case by using the shift operators which 
lower the multiplicity m a at one of the orbits only. 



Now let us return to Theorem 5.11 and apply the duality. Recall that starting 
from a root system R and multiplicities m a € Z + we have constructed a function 

We can cons ider also a similar function for the dual root system i? v . Then 
Theorem 4/7 implies that these two functions are related through interchanging the 
arguments: 

$-(ar, z) = $_(z, x) . 

Now take x = A S p m +i + P+ and z = /i£ Pm+i + ^+ an d use Theorem 5.11. This 
gives that 

c(A, m)^ m ( M )P A _ Pm+1 ( M ; g, <z" l+1 ) = c v ( M , m)^(A)P„ V __v ^ (A; g, g m+1 ) , 



-+1 



where we used V to distinguish the objects related to the dual root system P v . 
Taking into account formulas (5.6), ( 2.30| ) and their natural counterparts for the 
system i? v , we arrive after simple transformations to the relation 

P x . Pm+1 ( Ji ;q,q m + 1 ) J] JJ [(a v , A) + = 

P M _ p v +i (A;g,<z m+1 ) J] l[[(a,n)+j]-\ (5.10) 

aSi?+ ]=0 
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Put now p = pm+v then the right-hand side will not depend on A (since Pq(x) = 1). 
This gives us that 

Px- Pm+1 (pl + i,<l,Q m+1 ) = const I] l[[(a\X)+j}. 

aeR+ j=0 

To determine the constant, we substitute A = p m +i which leads to the evaluation 
identity: 

Denoting A — p m +i by v and m a + 1 by k a , we can rewrite it as follows: 

fc o _1 1 „27' + 2(Q V .^ + Pfc) 

p f/?/-<7 o fc ) - fl- 2 ^*'") TT TT g 



This can be presented as 



( a >")- 1 1 „2j+2fc Q +2(a v ,p fc ) 

^u{p k ,q,q )-q n n l _ 2j+2( a ^ Pk ) ■ 

aeR+ j=0 y 

An advantage of this form is that it works for all (e.g. non- integer) k a . Indeed, it 
is clearly rational function of t a = q ka . On the other hand, it is easy to see that 
for all k a the value P v (p^,;q,q k ) must be rational in t = q k . Since this identity is 
valid for integer k, it remains valid for all k. 

Using this identity (and its counterpart for the system P v ) one can rewrite the 
relation (5.10) in a more compact form, known as the symmetry identity: 



PA- Pm+1 G^+i; <?, Q m+1 ) P M V _ p v +i (p m +V,q, q m+1 ) 
Changing notations, we can rewrite it as follows: 
P x {p + pl-q,q k ) _ P^(X + Pk ;q,q k ) 



Px(pX;q,q k ) P^(p k ;q,q k ) 



VA e P+ , V/i e Pi 



Again, it is easy to see that both sides are rational functions of t = q k . Thus, in 
this form this identity is valid for all (e.g. non- integer) k a . 

This last identity is very important, since it leads directly (see [C2]) to the 
recurrence relations between Macdonald polynomials P\ with different A, which is 
a higher analogue of the three-term relation for classical orthogonal polynomials. 

5.6. Limiting case q = 1. Let us describe briefly what happens in the limit q, t — > 
1 if we keep t a to be of the form t a = q~ m * with fixed m a . It is convenient 
to present q as q = e e and rescale the rc-variable: x — > e x. The normalized 
BA function constructed in previous sections is a common eigenfunction of the 
Macdonald operators, acting in x and z. Let us rewrite these operators in new 
notations. 

First, in the a;- variable we have the operators D n related to (quasi) minuscule 
coweights of the root system R and given by formulas fl2.l|)~(f2!3|) with k a = — m a . 
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To avoid non-essential details, let us consider minuscule coweights only. Then in 
new notations the formula ( [2.1| ) takes the form: 

dz= e arir, o p( x)= n sinh( '.Tr +( r x)) - ( 5 - n > 

z — ' - LJ - sinh a.i 

tGWtt «EBi v ' ; 

(<»,t)>0 

Here T^ r denotes the shift in the x- variable by the vector er. 

Similarly, in the z-variable we have the operators related to minuscule weights 
lj of the system R: 

°"= j2 a ' T ' > n ""l^rr" - <« 2 > 

(0,T»0 

The normalized BA function ip(x, z) related to the root system R and multi- 
plicities m a £ Z+, is a common eigenfunction of all these operators. It can be 
characterized uniquely in terms of its analytic properties in the z-variable, as it 
was done above or, alternatively, in terms of its properties in x (due to duality). 
Similarly, one can construct it using the operator D" as in Theorem [T^] or, alter- 
natively, by using in a similar way the operator . 

Now let us look what happens in the limit e — ► 0. It is quite clear that the 
operator ( |5.12 ) in this limit takes the form: 

Z^=£a T TJ, a T (z)= J] ^J""" ' ^ 

(a , z) 



t£Wuj 

(o,t)>0 



On the other hand, it is slightly more difficult to see what is a proper limit of the 
operator (5.11) as e goes to zero. In order to do this one should expand in a 
series in e and pick up the first nontrivial term of this expansion. It turns out to 
be the following second-order differential operator: 

L m = A - 2m a coth(a, x)d a + 4(p m , p m ) , (5.14) 

where A is the Laplacian in V = M™, d a stands for the derivative in a-direction, 
and p m is given by (3.1). This operator plays the central role in Heckman-Opdam 
theory of multivariable hypergeometric functions [HO], and it is gauge-equivalent 
to the generalized Calogero-Sutherland operator from [OP|. 

From this it is natural to expect that a proper limit of the BA function ip(x, z) 
must give a common eigenfunction both for the rational Macdonald operators (5.13) 
and the operator L m . However, it is quite difficult to sec this directly, looking at our 
formula for if). The best way is to repeat the main constructions independently for 
this degenerate case. Since now we have no symmetry between x and z variables, 
we may choose two different ways to describe tp: either in terms of x-properties, or 
in terms of its properties in the z-variable. The x-part of the story is very similar 
to what we had before: ip has the form 

^{x,z) = e 2 ^Y,^ e2[U,X) 



with the summation taken over the polytope (4.21). The main difference is that the 
conditions (4.17) are replaced by the following: for each a £ R+ and s = 1, . . . , m a 



(<9 a ) 2s -V = for q 2 ^ = 1 

39 



(5.15) 



Here d a denotes the derivative in a-direction in the x-variable. One can check that 
the operator (5.14) with k a — —m a preserves these properties, which gives an ana- 
log of Proposition 2.1. After that everything becomes more or less a straightforward 
modification of the previous constructions. 

The z-properties of ip in this case change more significantly: it becomes quasi- 
polynomial in z in a standard sense: 

tp(x, z) = P(x, z) exp 2(x,z), P = C(x) J[ (a v , z) m " + ... , 

where C(x) = \\ a&R (sinh(ci!, x)) m " and the dots stand for lower terms in z. A 
detailed exposition of this latter approach is given in [Chi], where the relation of 
the ^-function to the hypergeometric function and Jacobi polynomials by Heckman 
and Opdam is also discussed in detail. 

To illustrate the difference between these two approaches, let us compare how 
the formula (3.16) will look for each of them. First, one can construct ip using 
L = L m . Let us introduce a trigonometric polynomial Q(x) as 

Q{x) = e 2 (^) Yl (2sinh(a,a;)) 2m ° , 



where p = p m is given by (3.1). Then ^i, up to a certain z-depending factor, is 
given by the formula 



4,^\{(L-A(z + uf)[Q{x)e 2 ^] 



where the product is taken over all v ^ having the form v — ^2 a>0 l a ct with 
l a 0,1,..., ttIq, . 

Alternatively, to construct ip one can use D = given by (5.13). Then, up to 
a certain z-depending factor, tjj is given by the expression (see Chi]) 



; (x)) M [q( Z )e 2 ^^}, M= m « 



where m ul (x) — X^revvw e 2 ^' 2 ^ and the polynomial q(z) looks as follows: 



aeR+j=l 



The first expression for tj) can be viewed as a trigonometric generalization of the 
original formula by Berest [Be], while the second one is its difference version. As 
we see, they both are specializations of the same formula ( 3.16| ). 



6. BC n CASE AND KOORNWINDER POLYNOMIALS 

In this section we consider the case of the non-reduced root system R = BC n . 
A proper generalization of the Macdonald theory in this case was proposed by 
Koornwinder |Ko2| . It depends on five parameters (apart from q) and generalizes 
Askey- Wilson polynomials |AW| to higher dimensions. The one-dimensional case 
will be essential for us, so we start considering R = BC\. 
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6.1. Rank-one case. In this case we have 4 parameters a,b,c,d (apart from q), 
and the corresponding one-dimensional difference operator D, suggested by Askey 
and Wilson, looks as follows [AW]: 



D = v+(z)(T- 1) +v-(z){T- 1 - 1) + q(abcd)- x + q^abcd , 



(6.1) 



where T* 1 denotes the shift by ±1 in z and the coefficients u ± are given by the 
following formula: 



v+(z) 



{aq z - a~ 1 q- z )(bq z - b~ 1 q- z )(cq z + c~ 1 q- z ){dq z + cT^ -2 ) 
(q z — q~ z )(q^ +z — q~ 



'■)(q z +q- z )(qi +z +q-2~ z ) 



v (z) =v + (— z) . 



(6.2) 
(6.3) 



Our notations differ from those of Askey and Wilson [AW]: what they denote by 
(q, a, b, c, d) is (q 2 , a 2 , b 2 , — c 2 , — d 2 ) in our notations. 

The following function A(z) = A(z; a, 6, c, d, q) plays an important role in Askey- 
Wilson's theory: 



A(z) = (abcd/q)~ 



(g 4^ 2+ 4 2;g 4 )c 



(a 2 q 2z , b 2 q 2z , — c 2 q 2z , —d 2 q 2z ; q 2 ) c 
Here we used the standard notations: 



(6.4) 



(*; «)c 



nt 1 -^ 



(ti, . . . , t„; q) c 



Using it, one can present the coefficients v of the Askey- Wilson operator D as 

v+(z) = A(z + 1)/A(z) , v~(z) = A(-z - l)/A(-z) . 
Introduce the dwaZ parameters a, 6, c, <i as follows: 

a=(abcd/q)^, b=(abq/cd)i, c—(acq/bd)^, d = (adq/bc)? . 
We will denote by A the function ( |6.4| ) with the dual parameters: 

A(z) = A(z; a, 6, c, d, q) . 



(6.5) 
(6.6) 



(6.7) 



Now let us make some special choice of parameters a, b, c, d in (6.2). Namely, we 



put 



— I j —l — m j 

a = q , o = g , c = g , d = q 



where l,l',m,m' £ iZ are some (half)intcgcrs with the requirement that 



^ + 1 + 1' £Z 



and — + m + m! £ Z . 



We will assume that l,l',m,m' are positive. Introduce iV as 

AT = l + Z + Z' + m + m' e Z . 



(6.8) 



(6.9) 



(6.10) 



Using it, one can present the constant term in (6J) as 



q{abcd) + q abed = q + q 
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-A? 



Below we will denote by (l,V ,fh,m!) the dual set of multiplicities, determined in 
accordance with (6.6): 



(6.11) 













i 


i 


i\ 


(i) 


V 




-1/2 


1 


1 


i 


-i 


-i 


fh 




-1/2 


+ 2 


1 


-i 


i 


-i 


m 










V 


-i 


-i 


V 


\m'J 



Notice that this transformation, as well as (6.6), is involutive. Geometrically, it 
reduces to the orthogonal reflection with respect to the hyperplanc l—l'—m—m' = 1 
in R 4 . 

Now we are going to formulate an analogue of Proposition 2.1 for the Askey- 
Wilson operator. First, introduce the following shorthand notation u =<! (v, w) in the 
situation when for 3 given numbers u,v,w at least one of the differences v — u,w — u 
belongs to Z> . Now let us consider a ring which consists of all polynomials 

jez 

satisfying the following N = 1 + 1 + 1' + m + m' conditions: 



f( z -s) = f(z + s) for 
f(z - s) = f(z + s) for 



;2Z 



q 2z = -1 



and 0<s4(lJ'), (6-12) 
and 0< a (m, to') ■ (6.13) 



Proposition 6.1. The Askey-Wilson operator (6.1)- fl6.3| ) with the parameters of 
the form (6.8)-(p^3) preserves the ring 9\ above: D(%K) C 9^. 



Proof can be found along the lines of our proof of Proposition 2.8. We need also 
the following inverse result. 



Let D be the Askey-Wilson operator (6.1) with a, 6, c, d as in (5.> 

2 cz where r is rational 



Lemma 6.2 

( 3.£ ). Let / be an eigenfunction of D of the form / 



r\x, z)q 



in z. Then / is, in fact, non-singular and satisfes conditions (6.12 )— (6.13) 



Proof. First, it is easy to see from equation Df = Xf that in the case ( |6.8| )-( |6.9|) 
any eigenfunction either has infinite number of poles, or has no poles at all. This 
proves that / is nonsigular in z. Then the conditions 
similar to lemma 2.7. 



.12 )— ( 5.13 ) can be derived 
□ 



All this suggests the following definition of a Baker-Akhiezer function for the 
Askey-Wilson operator with multiplicities (I, I', m, to'). 

Definition. A function tp(x, z) of the form 

V, = q 2xz Y, ^ - 1>i = 1>M , (6-14) 
b'l<iv 

is called a Baker-Akhiezer function for the Askey-Wilson operator if it satisfies 
N = l + l + l' + m + m' conditions (|T2|)-(|t|) in z for all x. 
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In the same way as in case R — A±, one proves that such a ip does exist and it 
is unique up to an x-depending factor. Also it is quite clear that the only nonzero 
coefficients ipj in ( |3.14| ) will be those with j + N G 2Z. However, we were not able 
to calculate ipj explicitly. This makes the following result somewhat less trivial. 

Proposition 6.3. Suppose that the multiplicities I, I 1 , m, mf € satisfy ( |6.9[ ) 

and are such that the dual multiplicities (6.11) are positive, too. Then if one puts 
tpN = A _1 (x) in accordance with the formulas (6.4)— (5.8), then all the coefficients 
ipj will be polynomial in x with ip—jsr(x) = ip]^(—x). 

As we already said, we cannot prove this proposition directly. So we use the 
following strategy: we will cons truct an cigenfunction ip of the Askey- Wilson op- 
erator which has the form (6.14) and satisfies the requirements of the proposition. 
The constructed ip will be automatically a BA function due to lemma p. 4 

To construct eigenfunctions of the operator (6T) for the case (3.8)-( |oT^ ), we will 
use shift operators, which in rank-one case can be computed directly. We need the 
operators which shift the parameters a,b,c,d in (6.1). We say that an operator S 
shifts from (a, 6, c, d) to (a, b, c, d) if the following intertwining relation holds: 

DoS = SoD, (6.15) 



where D and D are the Askey- Wilson operators (3.1) which correspond to (a, b, c, d) 
and (a,b,c,d), respectively. 

Proposition 6.4. Define A operators Si , S2 , S3 , S4 as 

S i =Af{z)T 1 '' 1 +Ar{z)T- 1 ' 2 , i = l,...,4, 
where Af (z) — A~l(— z) and Af are given by the formulas: 



At 


{aq z 


1 
2 


_i 1. 

— a g 2 


- z )(bq z - 


1 
2 


~b l q% 


~ z )(cq z 2 -f- c l qi z )(dq z 2 + d l q% z ) 












(q z - q~ 


- Z )(g z +<T Z ) 


At 


(aq z - 


1 

2 


_i I. 

— a q^ 


- z )(bq z - 


1 
2 


— b^ x q^ 








(q z - q- 


"*)(«* + 


(f 


~ z ) 




At 


{aq z ~ 


1 

2 


_i i- 
— a qz 


- z ){cq z ~ 


1 
2 


1 1 
■fc q 2 ~ 








(q z - q~ 


"*)(«* + 


T 


- z ) 


5 


At 


{aq z ~ 


1 

2 


_i i- 
— a L qz 


- z )(dq z ~ 


1 

2 


+ d~ 1 q^ 








(q z - q- 


- z )(q z + 


<r 


~ z ) 





Then 
shifts 

Si. 

S 2 - 

s 3 . 
s 4 



Si are the shift operators for the Askey-Wilson operator and the corresponding 
of the parameters are as follows: 

■ (a, 6, c, d) -> (q~^a, q~ib, q^^c, q'id) , 
(a,b,c,d) -> (q*a,q2b,q-2c,q-2d) , 

■ (a,b,c,d) -> (q2a,q-ib,qic,q~id) , 
(a,b,c,d) -> (qia,q-2b,q-2 C ,q2d) . 



Proof. It reduces to a straightforward though tedious calculation. Because of a 
certain symmetry between a,b,c,d it suffices to check the intertwining relation for 
Si and S*2 only. □ 



43 



Proof of Proposition \6.\ Recall that we assume that both 1,1' ,m,m' and their dual 



(3.11) are positive. Consider the case when I and fh are both integers. Then 
there is a proper composition of shifts Si which shifts from (1, q 1 ' 2 , 1, q 1 ' 2 ) to 
(q~ l , q~ 1 ' , q~ m , q~ m ')- Namely, one should apply Si I times, S2 {I' + 1/2) times, 
S3 fh times and S4 (fh' + 1/2) times. The Askey- Wilson operator with the pa- 
rameters (1, q 1 ! 2 , 1, q 1 / 2 ) is a trivial one, D Q = T + T -1 , with ^ = q 2xz being its 
eigenfunction. Applying the composition of the shifts above to tpo, we obtain an 
eigenfunction tji of the operator D. It is obviously quasipolynomial in x. More- 
over, it will be quasipolynomial in z, too (this is not completely obvious, since the 



shifts have singularities, but follows, for instance, from lemma 6.2). Now in order 
to find the leading coefficients, it is sufficient to look for each application of Si at 
the asymptotics of ip in z at ±00, using the explicit formulas for Si. After some 
simple inductive calculations, we arrive directly at the formula for i/)±n from the 
proposition. 

The three other possible cases, namely 
(l)I,m'6Z, (2)l',meZ, (3) !>'eZ, 

can be considered in a similar manner. The only difference is that in cases 2 
and 3 one applies shifts starting from D with (a, b, c, d) being (1, q 1 ! 2 , q 1 ^ 2 , q) or 
(1, q 1 ^ 2 , q, q 1 / 2 ), respectively. These operators are also almost trivial: they are 
obtained from Dq = T + T _1 by a simple gauge. For instance, the Askey- Wilson 
operator with parameters (1, q 1 ^ 2 , q 1 ' 2 , q) has the form 

D= q + q _ x (T-l)+ q +\ (T-l) + (q + q-l) 
qx _|_ g x q x + q x 

= Of + q-*)- 1 o(T + T- 1 ) o (q x + q- x ) . 

Then we apply Si (I + 1/2) times, S2 I' times, S3 fh times and S4 (fh' + 1/2) 
times, arriving at D with the parameters (q l , q l , q m , q rn ). Other arguments remain 
the same. The case 3 is analogous. □ 

This leads us directly to the main result of this section. 
Theorem 6.5. Let ^(x, z) denote a BA function with the parameters I, I', m, m' € 



Z + normalized as in Proposition 6.0. Then ^ solves the following bispectral system: 

D z ^> = (q 2x + q- 2x )^f, 



2: : , .2. ,,T, ( 6 ' 16 ) 



D x i/j = [q zz + q- Zz )$> . 



Here D z is the Askey-Wilson operator (6.1) related to 1,1' ,m,m! , and D x acts in 



x and is related to the dual multiplicities (6.11). Moreover, if the dual parameters 
are positive, we will have the duality as follows: *S>(z,x) = ^(x,z) where ty(x,z) 
denotes the normalized BA function related to the dual parameters. 



Proof. The first equation (in z) follows from Proposition 6.1 and the uniqueness of 
ip. It does not depend on our particular way of normalizing tp. To obtain the second 
equation we apply the standard argument: consider the function <fi = (q 2z + q~ 2z )^ . 



It is quasipolynomial in z, satisfying conditions (6.12)— (6.13). Then it must be 
obtained from ^ by applying a proper difference operator: 

L* = (q 2z + q- 2z )^ , L = b+(x)T x + MxjT- 1 + b (x) . 
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The coefficients b± can be found easily as soon as we know the leading coefficients 
v]/ ± jy in "J. To calculate 60, however, we need one more term, say, *&—n+2 in , 3 / . 
Unfortunately, we have not found anything better than to calculate it directly from 
the difference equation D z ^ = (q 2x + q~ 2x )"i> . This is pretty straightforward and 
we shall not reproduce this calculation here. As a result, one finds that the the 
operator L is nothing but the dual Askey- Wilson operator D x . 

Finally, t he duality between x and z follows similar to the case R = Ai, see 
Proposition 4.2. □ 



6.2. BA function for Koornwinder operator. Now let us consider the dif- 
ference operator by Koornwinder which generalizes the Askey- Wilson operator to 
higher dimensions. This oper ator D depends on five parameters a,b,c,d,t apart 
from q and it looks as follows [ Kol | : 

n j-2n y t~ 2n 

D = J2 v t(Ti ~ 1) +«rC*T 1 - + abcdq- 1 -— r + (abcd^q _ , (6.17) 

1=1 

where Tf stands for a shift by s in Zi and the coefficients Vf(zi, ■ ■ ■ , z n ) are given 
by the formulas: 

with the functions v (zi) obtained by substituting z — Zi into the formulas (5.2)- 
©■ 

The underlying geometrical structure here is an affine root system C v C n in 
notations of \M5\. For our purposes, however, it will be enough to consider a usual 
root system R of B„-type: 

R = {±e % } U {±e, ± e, \ i ^ 3} . 

We fix its positive half R+ as 

R+ = {ei, . . . , e„} U {ei ± e 3 \ i < 3} . 

The root lattice is Q = Z" and its positive part Q+ is defined as: 

3 

Q+ = {i/GZ"|^^>0, j = l,...,n}. 

i=i 

Weight lattice P in this case is also the standard lattice Z": 

P = {v = V\e\ H h v n e n \ Vi € Z} , 

while the cone of the dominant weights P+ looks as follows: 

P+ = {v e P I Vi > v 2 > ■ ■ ■ > v n > 0} . 

The Weyl group W acts by permuting the variables and flipping their signs arbi- 
trarily. 

As before, by a polynomial we mean a finite sum f(z) = J2, yeP q 2 ^ 2 ) , keep- 
ing calling functions like f(z)q 2 ^ x,z ^ quasipolynomial in z. Algebra of VF-invariant 
polynomials in this case is a linear span of orbitsums 

m x (z) = J2 9 2(T ' 2) . ( 6 - 19 ) 

t£W\ 
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2z, 



-2zi 



and it is a symmetric polynomial algebra of the generators = q 
Now let us specialize the parameters a, 6, c, d, t as follows: 

t = g~ fc , o = g _I , b = q- 1 ', c = q- m , d = q~ m ' , (6.20) 

where k G Z + and Z, Z', m, m! £ are some (half) integers with the requirement 
as in rank-one case that 



- + 1 + 1' €Z 



and 



— + m + m € Z . 



(6.21) 



We will denote by M the whole set 

M = (k,l,l',m,m') . 
Below we will also use the multiplicities m a defined in the following way: 



k for a = ±ei ± e,- 



and 



to q = l + Z + Z' + m + m' for a = ±e, . 
Introduce a vector p depending on M = [k, I, I', to, to') as follows: 



P= PM 



2 ^ 



(6.22) 
(6.23) 

(6.24) 



Using it, one can present the constant term in ( 6.17| ) in the following way: 

7 2(p,±e0 



abcdq^ 1 — — + (abcd)~ 1 q- 



1-f 



1 -t~ 



5> 



which makes it similar to the constant term in (2.2). 

Let 91 denote a ring which consists of all polynomials f(z) — X^es 
the following properties: 

(1) for alH = 1, . . . , n and < s =4 (I, I') 

(T? ~ Tr s )f = for q 2z ' = 1 ; 

(2) for alH = 1, . . . , n and < s =^ (m, to') 

(T° - Tr»)f = for^=-l; 

(3) for all 1 < i < j < n and s — 1, . . . , k 

(T? - T/)/ = for q 2z > = q 2z > ; 

(4) for all 1 < i < j < n and s — 1, . . . , k 

{T? - Tr s )f = for q 2z * = q~ 2 ^ . 



f(^-) with 

(6.25) 
(6.26) 
(6.27) 
(6.28) 



We used the same notation s =^ (Z,Z') as above, denoting that at least one of the 
differences I — s, V — s is a nonnegative integer. 



Proposition 6.6. The Koornwinder operator ( |6.17 )-( 6.18| ) with the parameters 
( 3.2C )— (6.21 ) preserves the ring 9t as above: D(9l) C 91. 

Moreover, in a similar way one can check that all n commuting difference opera- 
tors Di, . . . , D n constructed in |vDl], will preserve the ring 9t in case ( 6.20] )-( |6.21 ). 

Now a Baker- Akhiezer function ip(x, z) is defined similarly to the case of a 
reduced root system: 
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(1) ip has the form 
where the summation is taken over v G P lying inside the polytope 



N = {u= Li 



1 1 , 

'2 m a < L < 2 m a/ I 



(6.29) 



(6.30) 



+ 



(2) ip(x,z) satisfies the conditions (6.25)— (3.28) in z for all x. 

One proves, similar to Section |[ that such ip does exist and is uniqu e up to an 
x-depending factor. It can be expressed by a formula, similar to (3.16). Namely, 
introduce first a function c(A, M ) depending on A 6 R™ and the multiplicities 
M = (k, I, I', to, rn') in the following way: 

C (a,m) = n ii - ( a > a )] n b'-K A )i n -(«><*)]. ( 6 - ai ) 

QEijO n£H| o£S| 

0<j^(i,i') 0<j=$(m,m') j = l,...,fc 

where R and i? 1 denote the sets of short and long roots, respectively, and [a] 
stands, as usual, for q a — q~ a . 

Now define the polynomial Q(z) as follows: 

Q(z) = (-l) n i l + l '+^q 2 M c (z, M)c(-z, M) , (6.32) 

where c(z, M) is defined above and the vector p is given by ( 6.24 ). Introduce also 
the notation m for the orbitsum 



m(x)=^(^+g- 2 ^)- 



(6.33) 



Theorem 6.7. Let D be the Koornwinder operator (6.17) with the parameters as 
in ( 3.2C )-( 3.21 ) . Define ipix, z) as follows: 



J2{x,z) 



(6.34) 



in accordance with the formulas ( 3.32] ), ( |6.33 ), where the product is taken over all 



v 7^ having the form v = ~^2 a€R+ l a o^ with l a = 0, . . . , m a . Then 

(i) ip is a BA function for the Koornwinder operator; 

(ii) the coefficient ip— p in its expansion ( 6.29| ) equals 



H(m(x)-m(x + v)) ^0; 

V 

(Hi) as a function of z, ip is an eigenfunction of the Koornwinder operator D: 
Dip = m(x)ip. 



To normalize ip, consider the dual parameters 

M = (fc, I, l\ rh, to ) , 



where k = k while other 4 parameters transform according to (6.11). Let us nor- 
malize ip in the following way: 



ip p (x) =c{x,M) 
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(6.35) 



in accordance with the formula (3.31). 



Theorem 6.8 (Duality). Let ip(x,z) be a BA function related to the parameters 
I, V, to, to', k and normalized as above. Suppose that the dual parameters I, I', m,m',k 
are positive. Then tp(x, z) is quasipolynomial in both x and z and has the following 
duality property: 

ip(z,x) = ip(x,z) , 

where ip denotes the normalized BA function related to the dual parameters. In 
particular, ip(x, z), as a function of z, is an eigenfunction of the Koornwinder 
operator ( 6.17| ), while in x it satisfies a similar difference equation related to the 
dual parameters. 



The algebraic integrability of the Koornwinder operator in case (6.20) and the 
existence of the shift operators is a straightforward generalization of the similar 
results for reduced root systems. 

6.3. Koornwinder polynomials. The Koornwinder polynomials 

P\(z) = P\(z;q,a,b,c,d,t) 

can be defined similar to Macdonald ones, as polynomial eigenfunctions of the 
Koornwinder operator ( |6.17 ), see [Ko2. vDl]. Since our notations are slightly dif- 
ferent, we reproduce here their definition for the reader's convenience. Namely, P\ 
has the form 

P\ = vcix + ^2 ax v xa v , A G P+ , (6.36) 

in notations of the previous section, with v -< A meaning that A — v G Q + . For 
generic values of the parameters a, b, c, d, t the polynomial Pa is uniquely determined 
from the equation: 

DP X - c AA P A , (6.37) 
where D is the Koornwinder operator ( |6.17| ) and the eigenvalue c\\ has the form 

~ 7 2 (^+p), (6.38) 



with p = pm given by formulas ( 6.20| ) and (6.22 )— ( 6.24 ) 



To compare with the notations in [vDl|, one should put q = e lu} , then our 
orbitsums ttia(^) correspond to m\(ioz) in notations of [vDl|. Our defining equation 
(3.37) corresponds to the equation 



Dip\ = Ei >n {\ + p)p\ 



which is the case r — 1 of eq. (3.72) from [vDl] 



Now the relation of ip to Koornwinder polynomials (generalized Weyl formula) , 
the norm formula, evaluation identity and duality, they all can be derived similar 
to the case of a reduced root system. To avoid repetitions, we skip the details (see 
|vD2| ] for the formulation of all these identities) . 
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7. Integrable deformation of the Macdonald-Ruijsenaars operators 



In this section we discuss a version of the Macdonald-Ruijsenaars operators 
related to a certain "deformed" A n system. This system R — A n (m) C R n+1 was 
introduced in [CFV1] in the following way: R — R° U R 1 where 

R° ={±(e, - ej) | 1 < i < j < n} S A n _i C R" , 

i? 1 ={±(ej — y/me n +i) | i = 1, . . . , n} with (7-1) 

m a =m for a S R° , m Q = 1 for a G i? 1 . 

Here m is a parameter, which at first will be an integer. 

7.1. Deformed Macdonald Ruijsenaars operator. Let us consider the follow- 
ing difference operator related to the deformed A n system: 

D = a 1 (z)T ei + ■ ■ - + a n {z)T e " + a n +i(z)T^ e " +1 , (7.2) 
, e„+i is the standard basis in V = R n+1 and the coefficients look 



where e\, 
as follows 



di(z) 



a n +i(z) 



mz n+ i 



m+l 1 
2 J 



mz n+ i\ 



[ml 11 



n 

3^1 



mz n+1 - Zj 



[Zi 



m 



(7.3) 



m+l l 
2 J 



mz n+ i 



2 J 



Here, as before, the square brackets are used to denote [a] := q a — q~ a . 

This o perato r is a discretization of the deformed Calogero-Moser operator pro- 
posed in |CFV1|, with m being the deformation parameter. Its rational version was 
considered in [Chi, |Ch2 |. Fo r m — 1 it reduces to a special case of the Macdonald- 
Ruijsenaars operator D\ in (2.7). 



We fix a "positive half of the system (7.1) as 

R+ = {e t — ej | 1 < i < j < n} U {e« - Vme n+ i | i = 1, . . . , n} . 

Next, we introduce "weight lattice": 

P = Zei © • • • 8 1e n © -sfmlen+x . 

Respectively, we will call any finite sum f(x) = X)i/eP fvQ 2 ^' x ^ a polynomial in x. 
We don't need coroots and coweights in this case, and there will be no substantial 
difference between variables x, z e M ,l+1 below, so, for instance, polynomials in z 
arc defined in the same way. As before, we will apply the term quasipolynomial to 
a function ?p(x, z) — q 2 ( x ' z *) f where / is polynomial either in x or in z. 

Now let us consider a ring £H of polynomials f(z) with the properties ( |3.9| ) (with 
a v := a). In our case they can be rewritten as follows: 

f(z + se % ) = f(z + se j ) for q 2zi = q 2z > (1 < s < m, 1 < % < j < n) , (7.4) 



f(z + et) = f(z + V^e n+1 ) for q ^ m +^ = gV^+i (j = l, 



(7.5) 



Proposition 7.1. Deformed Macdonald-Ruijsenaars operator ( 7.2 ) — ( |7.3| ) withm € 
Z + preserves the ring 91 o/ f/ie polynomials with the properties ( 7.4 ) — ( |775| ) / _D(£H) C 

4!) 



Proof. For the conditions J7.4| ) the arguments repeat those from the proof of Propo- 
sition 2.1, because the operator D is symmetric in z\, ■ ■ . , z n . For the remaining 
conditions (7.5) everything reduces to the proof that the operators Dq = a,j(z)T ej 
(for j ^ i) and D 1 = a,i{z)T e ' +a n+1 (z)T^ Ke ^+ 1 preserve the property fl7l| . This 
can be checked straightforwardly, in the spirit of Lemma 2.2. □ 



7.2. BA function. Introduce p — p(m) similar to (3.1): 

P = - ^2 ma01 - 



(7.6) 



u£R 



+ 



In our case we have explicitly: 
m 

P= y( n ~ l,n-3,.. 
We define a polytope N similar to 

N = {u = -p- 



-n + l,0) + -(l,. 



ran) . 



l a a\0<l < m a } ■ 



(7.7) 



Now the definition of a BA function repeats our definition in case of a root system. 
Definition. A function ip{ x : z ) which is quasipolynomial in z, 



2{y,z) 



ueP 



with supp(^) C J\f and which satisfies the conditions (7.4)- ( |7T5| ) in z is called a BA 
function for the system (7.1). 

In exactly the same way as in Section |§|, one proves that tp is defined uniquely up 
to an x-depending factor. Moreover, analyzing the corresponding linear conditions 
for the coefficients ip„, we come to the following choice of normalizing tp: 

^= n ii^2 ja ~ x)h (7 - 8) 

with [a] denoting the g-numbcr [a] := q a — q~ a . We will call this ip the normalized 
BA function. One can calculate the leading coefficients ip v at other vertices of the 
polytope jV. It has (n+l)l vertices which are in one-to-one correspondence with all 
permutations a £ S n +i- Namely, for a € S n +i let us introduce a vector t a £ M. n+1 
as 

t a = (<7i, . . . , a n ,cr n+ i/y/rn) . 
Now denote by aR + the following subset in R: 



and introduce op as 



aR + = {a £ R | (a, t a ) < 0} , 
- ^ m a a. 



op 



a£t?R 



+ 



These are exactly the vertices of the polytope (7.7). For instance, taking o — ewe 
get op = p, and taking o = (n + 1, n, . . . , 1) we get op = —p. 



Similarly to Proposition 4.4. one gets the following result. 
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Proposition 7.2. The leading coefficients ipo-p of the normalized BA function have 
the form: 



irt/ot i 



(7.9) 



The normalized BA function is quasipolynomial in both x and z. 



The existence of ip can be proven similar to Theorem 3.7. To formulate the 
result, introduce a polynomial Q(z) as 



1 1 

Q(z) = q 2 ^ J] l[[(a,z+-ja)}[(a,z--ja)}, (7.10) 

a£R + j=l 

with square brackets denoting g-number as before. Define also a "deformed orbit- 
sum" m(x) as follows: 



m(x) = q 



2.ri 



q m - q- 



2^/rriXr, 



(7.11) 



Theorem 7.3. Let D be the deformed Macdonald-Ruijsenaars operator ( [7.2| )- 
([7.3|). Define ip(x, z) as follows: 



ip = Y\(d - m(x + v) 



,2(x,z) 



Q{z) 



(7.12) 



in accordance with the formulas ( 7.10] ) ; (7.11), where the product is taken over all 
v =/= Q having the form v — ^2 aeR l a ot with l a = 0, . . . , m a . We have the following: 



(i) %j) is a BA function for the system (7.1); 

(ii) its coefficient ip- p equals Y[u( m ( x ) ~ m ( x + v )) 7^ 0/ 

(Hi) as a function of z, ip is an eigenfunction of the operator D: Dtjj = m(x)ip. 

Thus, renormalizing the constructed if) one gets the normalized Baker- Akhiczcr 
function \& for the system (7.1). 



Now one can derive the duality similar to Theorem 4.7. In this case it is simply 
the symmetry between x, z. 

Theorem 7.4. The normalized BA function ^ constructed above, is symmetric 
under permutation of its arguments: ^(x,z) — ^>(z,x). 

7.3. Quantum integrability. We start from discussing the algebraic integrability 
of the deformed Macdonald-Ruijsenaars operator ([7.2[)— (ffT^) with m e Z + , which 



is a direct corollary of the existence of a BA function for the system (7.1 



Theorem 7.5 (Algebraic integrability). Let m 6 Z + and ^{x, z) be the normalized 
B A function for the deformed A n system (7.1). Then for each polynomial f (z) from 
the ring 91 there exists a difference operator Df in x on the "weight" lattice P such 
that Df'fy = f{z)^. All these operators commute. For f = m(z) € SH given 
by (7.11), the corresponding operator Df is the deformed Macdonald-Ruijsenaars 



operator D given ( p.l| ) -(2.3) 



Proof is the same as in Theorem 5.1. In this way we obtain a commutative ring 



of difference operators, isomorphic to the ring (7.4)— (7.5). Moreover, due to the 



symmetry between x and z, we obtain a bispectral pair of commutative rings with 
^{x,z) being their common eigenfunction. 

The ring is big enough, it contains, for instance, a principal ideal gener- 



ated by the polynomial ( 7.10 ). Now let us choose the following special elements 
mi, ... , m„ + i of the ring d\: 



s 



m s (z) = q 2sz > +■■■ + q 2s *" + J- 9 q 2 ^ z " +1 . (7.13) 



In particular, for s — 1 we obtain the deformed orbitsum from (7.11). Notice that 
for m = 1 these polynomials turn into the Newton basis in the ring of symmetric 
functions. 



In accordance with Theorem 7.5, to each m s corresponds a certain difference 



operator D s = D mg , and they all commute. Similar to Lemma 5.5, one can show 
that these operators D s admit analytic continuation in m and, thus, they give rise 
to a commutative family for any value of the parameter m. Thus, we arrive at 
quantum intcgrability of the deformed Macdonald-Ruijsenaars operator. 

Theorem 7.6. The operator D is completely integrable, i.e. it can be included into 
a commutative family D\ — D, D2, ■ ■ ■ , -D n +i of difference operators, which in case 
m = 1 coincide with the Macdonald-Ruijsenaars family ((2.7|). 



One obtains a natural elliptic version of the operator (7^2) replacing all expres- 
sions like [a] = q a — q~ a in its coefficients by their elliptic analogues a(a) with 
<j(z) = <j(z\bj,bj') being the Weierstrass c-function. For such an operator a proper 
version of Proposition holds. This indicates that it is (algebraically) integrable, 
too. As a concluding remark, we mention that in a similar manner one can con- 
struct generalized Macdonald operators for other "deformed" root systems, some 
of which were presented in JCFV2 |. Details will appear elsewhere. 
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